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TWO- AND THREE-DIMENSIONAL UNSTEADY LIFT PROBLEMS 
IN HIGH-SPEED FLIGHT ' 


By Harvarp Lomax, Max. A. HEABLET, FRANKLYN B. FULLER, aud Loma SLUVER 


SUMMARY 


The problem of transient lift on two- and three-dimensional 
wings flying at high speeds is discussed as a boundary-calue 
problem for the classical ware equation. .Kirchhoff's formula 
is applied so that the analysis is reduced, just as in the steady 
state, to an investigation of sources and doublets. The appli- 
cations include the evaluation of indicial lift and pitching- 
moment curres for two-dimensional sinking and pitching 
wings flying at Mach numbers equal to 0, 0.8, 1.0, 1.2, and 
2.0. Results for the sinking case are also given for a Mach 
number of 0.5. In addition, the indicial functions for super- 
»onic-edged triangular wings in both forward and reversed flow 
are presented and compared with the two-dimensional values. 


INTRODUCTION 


The usual idealizations introduced in the development of 
linearized aerodynamic theory describe a frictionless, per- 
feetly elastic, model fluid. As is well known, the effect 
of small disturbances in such a fluid can be analyzed by 
means of the familiar wave equation. 

The solution to the wave equation known 8s Kirchhoff's 
formula (see reference 1) is found to be of considerable use in 
unsteady-motion problems involving thin wings with super- 
sonic edges. The problem is reduced to one of summing 
elementary solutions, analogous to sources and doublets 
in steady flow, over a region determined by the position of the 
Wing as well as its traversed path. The theoretical develop- 
ment leads naturally to the concepts (defined later) of, first, 
inverse sound waves, whieh have a counterpart in the Mach 
forecones used in steady-state wing theory; second, acoustic 
plan forms; and, third, homogeneous flow, which reduces 
in part to the familiar conical flow as the wing approaches 
a steady supersonic velocity. 

There are several simple types of unsteady motion on 
Which the analysis can be based. The so-called indicial 
motion, in which the velocity undergoes an initial discon- 
tinuous change will be considered here. It is possible to 
conceive the perturbation field due to the unsteady motion 
in two slightly different ways. For one, it can be supposed 
that the wing has been traveling at the constant velocity 
V, for an infinitely long time and then, at time equals zero, 
starts suddenly to sink without pitching (or pitch at a 


constant angle of attack) while maintaining the forward 
velocity Wy. On the other hand, the wing may be considered 
to be at rest in still air until at time zero it starts suddenly 
either to sink or to pitch and, at the same instant, attains 
the forward velocity W>. The latter physical picture will be 
used in this report. Problems of unsteady motion can also 
be approached with the initial assumption that the velocity 
potential depends harmonically on the time. These two 
approaches are quite compatible in that they can be related 
through the use of superposition methods (Duhemel’s 
integral, Fourier’s integral) of the operational calculus. 
Detailed results for two classes of indicial responses (time 
response to a step input) will be given in this report. The 
first of these is the complete set of responses (ei. Cmar Cte: Eme) 
for two-dimensional wings flying at Mach numbers equal 
to 0, 0.8, 1.0, 1.2, and 2.0. In addition, results for e: and 


Cm, are given for a Mach number equal to 0.5. The part of 


the analysis pertaining to the response of wings traveling at 
subsonic speeds is lengthy and somewhat tedious regardless 
of the method of approach. With the use of indicial func- 
tions, however, the calculations are reasonably straightfor- 
ward, especially for Mach numbers around 0.8 to 1.0. 
Further, the use of indicial functions sheds considerable 
light on the manner in which Mach number variations 
affect the section aerodynamic characteristics. The results 
for the two-dimensional wings traveling at supersonic speeds 
are compared with the indicial responses developed by 
rectangular wings of aspect ratios 2, 4, and 6, the latter 
curves having been presented by J. W. Miles in reference 2. 

The second class of indicial responses considered is that 
concerned with the forces and moments induced on sinking 
and pitching triangular wings. First, the loading on a flat 
triangular wing with supersonic edges undergoing an indicial 
sinking motion is determined. (See also reference 3.) Then 
a simplified method is developed whereby total lift and 
pitching-moment coefficients for the wing with supersonic 
edges may be obtained. (See also reference 4.) Again a 
complete set of indicial responses is presented for supersonic- 
edged wings in both forward and reversed flow. Lasily, 
the triangular wing with subsonic edges is partially analyzed, 
and the indicial responses for a slender triangular wing are 
given. 


| Supersedes NACA TN 2403, “The Indicia! Lift and Pitehing Moment for a Sinking or Pitching Two-Dimensfonsl Wing Flying at Subsonic or Supersonic Speeds” by Harvard Lomax 
Mar. A. Heaslet, and Loma Sluder, 1951; NACA TN 2357, “Three-Dimensional Unsteady Lift Problems in High-Spead Filght—The Triangular Wing” by Harvard Lomar, Max. 4. Heaslet, 
and Franklyn B. Fuller, 1951; and alse contains material from NACA TN 2256, ““Three-Dimensional, Unsteads-Lift Problems in High-Speed Flight—Basie Concepts" by Harvard Lomax, ` 


Max. A. Heasiet, and Franklyn B. Fuller, 1951. 
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SYMBOLS 


speed of sound in the undisturbed fluid 
chord of wing 


lift coefficient f 4- an 
z Po Vos 
indicial lift coefficient due to angle-of-attack 


change, without pitching (o. RE : ) 
indicial lift coefficient due to pitching for a 
wing rotating about its leading edge or apex 


Vo ) 


pitching-moment coefficient, positive when trail- 
ing edge tends to sink relative to leading edge 


moment 
] - 
5 Po V ges 


indicial pitching-moment coefficient due to angle- 
of-attack change (without pitching) measured 
about the leading edge or apex; positive when 
trailing edge is forced downward with respect 


to the leading edge or apex (Cn = ) 
indicial pitching-moment coefficient due to pitch- 
ing measured about the leading edge or apex, 
for a wing rotating about its leading edge or 
apex; positive when trailing edge is forced 
downward with respect to the leading edge or 


apex | C, = OC. ] 
"|a. 


two-dimensional lift coefficient e 
2 po Voc 

















two-dimensional pitching-moment coefficient 


moment 

(: Po ve) 
cotangent of sweep angle (cot A) 
Mach number in the undisturbed fluid 
M. rans 
2a8 Jus Yo y 
m Mo Vo Í "Apa 

268? —s Qo 
loading coefficient (pressure on the lower surface 

minus pressure on the upper surface divided by 


free-stream dynamic pressure) 


1 
free-stream dynamic pressure (z poVit) 


dimensionless rate of pitching (+ 
0 

local semispan of wing 

Wing area, 

time 

at^, aot’ fe 


u,v, W perturbation velocities m z, y, 2 directions, re- 
spectively 
V, . free-cstream velocity 
t, Y, 2 Cartesian coordinates, fixed relative to the fluid 
at infinity 
(z/c)..;. distance of center of pressure from leading edge 
or apex, in percent chord 
a . angle of attack (angle between flight path and 
plane of wing), radians 
| do 
» di 
B vii — 4] 
r cireulation 
S(t’) Dirac 6 function, normalized with respect to t, 
thus à(£/)-—0 for ¢’+0, 8(0) — o, and 
[| sari 
A discontinuity in the quantity in question across 
the plane of a wing 
8 wing angle of pitch relative to horizontal, posi- 
tive when trailing edge lies below leading edge, 
radians 
ZN d 
wing rate of pitch dp 
A angle of sweep of leading edge, positive for sweep- 
^ A back 
Po free-stream density 
: 
To -© chord lengths traveled (= ) 
e perturbation velocity potential 
| SUBSCRIPTS 
u upper surface of a wing 
i lower surface of a wing 
z,y,2,€ indicate differentiation with respect to the vari- 


able in question 


PART I—THE USE OF INDICIAL FUNCTIONS 
IN UNSTEADY LIFT PROBLEMS 


In the first part of this report, a discussion of three un- 
steady lift problems will be given. The three problems are: 

(i) Determination of the indicial response in lift and 
moment on a sinking or rotating wing; 

(ii) Determination of frequency response in lift and mo- 
ment on & fluttering wing; and 

(Àj) Determination of frequency response in lift and mo- 
ment on & slowly oscillating wing. 


The method by which the latter two problems may be 
solved with the aid of the solution to the first is also de- 
scribed, as well as the application of these results to the de- 
termination of lift and moment on wings undergoing arbitrary 


maneuvers. 
THE INDICIAL FUNCTIONS 


By definition, an indicial function is the response to a 
disturbance which is applied abruptly at time zero and is 
held constant thereafter; that is, a disturbance given by a 
step function. For example, if the angle of attack of a wing 
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varies with time as shown in figure 1, the corresponding re- 
sponse, also shown in figure 1, is designated as the indicial 
lift coefficient due to angle of attack. Four such indicial func- 
tions will be evaluated, namely, Cra, Cx,/, Cu’, Cu’. The 
primes on the coefficients indicate that the axes about which 
pitching motion occurs and pitching moments are measured 
either coincide with the leading edge of the wing or with a 
line through the apex normal to the root chord of the wing. 

The equations which transform these functions to those 
for a wing pitching about an axis a distance ac back from the 
leading edge and having its moment center a distance be back 
from the leading edge are simply 


Ce Z Cu + b Cr: 
CL = Cy a Cz 
Um m n + b CL'—a Cra —ab Cia 


a 


JT 





Fiecgz t.—Indicial excitation and response. 


In linearized, thin-airfoil theory the boundary condition 
that applies to the indicial functions due to angle of attack 
is simply 

Wa =W= — Vow (2) 


over a certain planar area in the ryz space at the time /'. 
If we consider & coordinate system fixed relative to the fluid 
at infinity and a wing with span measured along the y axis, 
moving away from the origin along the negative x axis with 
a velocity Fo, then equation (2) applies to the area in the 
z= plane occupied by the wing at any given time. It is 
further required that the ¢ be continuous everywhere except 
across the wing plan form and wake. In the case of the 
indicial functions due to pitch for a wing rotating about its 
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leading edge or apex, the boundary conditions for a wing the 
leading edge or apex of which is at the origin at /^—0 are 
that the upwash be given by the expression | 


1p, 1p; — (z+ Vol’) 6 (3) 


over the same region in the z—0 plane as for the angle-of- 
attack case and, again, that » be continuous except across 
the wing and its wake. The angle of pitch, 0, is taken as 
positive when the trailing edge is lower than the leading edge, 
and Å is the time derivative, d6/di', positive when the trailing 
edge is falling with reference to the leading edge. 

The difference between 8 and a is illustrated in figure 2 (a). 
The angle of attack o is the angle between the flat wing 
surface and the tangent to the flight path of some point 
fixed at a distance ac back from the wing leading edge. For 
example, in applications involving the dynamic behavior of 
an entire airplane the distance ac would usually be taken as 
the distance back to the center of gravity of the airplane. 
The angle 8 is the angle between the Hat wing surface and 
the horizontal. Figure 2 (b) shows a wing undergoing a 
sinusoidal angie-of-attack variation with & zero angle of 


Direction of wing mofion 









Fiigh? path of reference 
point on wing--~ 


bom 
sn 
(8) 
ge E ra 
— <= : 
iU NS 8-0«-0 
— M e-sin wt 
> x 
Sas 
— a — a 
———— Pee 27 
(b) 
A e ii 
Padi E 85sin wt 
PEL 
M 
c 
C "di 

(c) 

(a) Angle convention. 

(b) Sinking wing. 
(c) Pitching wing. 


FIGURE 2—Definition of angles of attack and pitch. 
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pitch throughout. For convenience and in order to dis- 
tinguish from the pitching wing, such a wing will be referred 
to as a sinking wing. Figure 2 (c) shows a wing undergoing 
a sinusoidal angle-of-pitch variation taken about the axis ac 
back from the leading edge. In the case shown the wing has 
a constant (zero) angle of attack. - 

The variable g to which the lift and pitching-moment 
coefficients of the pitching wing are referred is equal to 
éc/Vo, & dimensionless form of expression for rate of pitch. 


ADAPTATION OF THE INDICIAL FUNCTIONS TO MANEUVERS 


Since the theory is linear, the lift and moment on & wing 
undergoing an arbitrary maneuver Gn which both œ and 6 
remain small) can be calculated from the jndicial functions 
by the principle of superposition. For example, consider a 
wing performing some arbitrary maneuver involving small 
variations of both æ and 6. Let a and @ be defined ? with 
respect to the flight path of a point ac back from the leading 
edge. The appropriate indicial functions can then be cal- 
culated from those referred to the flight path of the leading 
edge by means of equations (1). (The moment center is 
still arbitrary at & point a distance be back from the leading 
edge. However, for convenience, b is usually taken to be 
equal to a.) Finally, by Duhamel’s integral, the lift and 
moment induced on the wing by the arbitrary maneuver are 
given by the equations 


Ef 
C. ds (1C, CHa) + On, UH) 1) at 
dc (4a) 
Om (10, 80 a+ On, P600 62) dt 
or alternatively 
d 2 / / f 
Cog f [Cr at iH O tE —5]dts 
(4b) 


du a | 
Cem | Ico (eral H+ Ce ge bah 


where, for example, C ea! ) meens that the indicial 
function Cy, is to be evaluated at the time t —&'. 


THE FLUTTER DERIVATIVES 


By definition, & flutter derivative is the response to & 
harmonie oscillation in vertical displacement (angle of attack) 
or pitch, such oscillation having continued so long that all 
transient effects of its origin are damped out and the induced 
forces and moments are periodic. In the notation used in 
reference 5 the induced responses are given in terms of the 
flutter derivatives Li, La, Ls’, Le and Mt’, My’, Mi, AL. 
The relation between the indicial responses and the flutter 
derivatives follows from equations (4) and can be written 
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: iett d ft 
Wie ne E I P an (Eiet dt 
i'-o 2k t 0 "a * ; 
L,-àiL/-—lim Rond a E [c (t^ —t ^4- 
t'—o Ak? di 0 ta 3 (5) 


2ike, (t —t)]e*^'dt,' 


—iíot 


pai r . —É€ d ct f ; 
Af,’ tiM: = lim a) ae], [Cm EC —45^) + 
2ikes ' (t —6)]e ^^ dt 


where e is the frequency and & is the reduced frequency, 
we/2V. The primes on the quantities again indicate that 
the wing is pitching about and the moments are measured 
about the leading edge of the wing. 

The responses in lift and moment. for harmonically oscil- 
lating wings can be expressed in terms of their absolute 
values and phase shifts with reference to the period of the 
forcing disturbance. Thus, if an angle of attack variation 
given by the equation o—e'*" is impressed on the wing, the 
response in lift would be le,|é“"*". In terms of the 


indicial functions this can be written 


hi f d e f 
le; lef* — lime 75] ei (t — thet dt,’ 
t'—o 2 


dt’ Jo 


which, by comparison with equations (5), can be re-expressed 
in terms of the flutter derivatives as 


Jer Jett =4k(L,—1L,) = 4k JLF Lye tn H) 


Similarly, each of the stability derivatives can be expressed 
in terms of the indicial functions and the flutter derivatives, 
hence: . 





Term Maximum value Phase angle 
€i, Ak L4 TLLA —ian^ 7: 
- 
" t 
Cm” 2k AM P+M tan"! My 
z Ma 
| LIO 
ik(ei;-- 0i) 2k4 L - L4 tan”! Ly 
ik (emg Hem) E M? M tan" Am 
$ 


where it is assumed that unit absolute values of angle of 


i : ; . _dajdť 
attack and pitch are impressed and the term à-c 777 —* 
7 Ü 


+ Notice that if a' and 6" are the angle of attack and pitch measured with respect to the flight path of the leading edge, and a and 6 the same angles measured with respect to the flight path 


ofa point ac back from the leading edge, the relation between the two seta of angles is given, 


for small deflections, by 


a! =a 
vi 
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A WING UNDERGOING SLOW HARMONIC OSCILLATIONS 


The boundary conditions for these problems are exactly 
the same as those used for evaluating the flutter derivatives, 
being in fact a special case of the latter when the frequency 
of oscillation is very low. This added simplification is im- 
portant, however, because the frequency of angle of attack 
and piteh oscillations for the entire airplane is, relative to 
the flutter frequencies, very low. Further, special methods 
based on the assumption that the wing is oscillating slowly 
have been devised. Thus, if the lift and moment are ex- 
panded in powers of the frequency, solutions for the coef- 
cients of the lowest-order terms have been found for both 
triangular and rectangular wings-flying at supersonic speeds. 
(See, e. g., references 6, 7, 8, and 9.) 

For a given set of indicial functions these expansions can 
he carried out with relative ease. Consider the case of Li 
and L as given bv equations (5). If the indicial curve for 
e, (t’) contains no pulse function,’ the expression for Lii, 
can be written 


Int iL lin ' — aN [ei (e ce )— Ac; ,U)]e gie -in di! 


where Ac; Q^) is the difference between the indicial Hft and 
its asymptotic value, thus 


Ae, UW) =e1,( 0) — ei P) 
Since Ae; (@) is zero, it follows 
Loi E (2) ia f Ac; Geze de | 
4 f a 0 a 


or, in terms of the reduced frequency parameter k=ac/2V4 
und the number of chord lengths traveled ra= Vt’ /e, 


L4 iL Cial »)-2ik| xe raed | 


In esse the flow is supersonic, these expressions can be 
expanded in powers of &, thus 


1L,—2 |, Aciatreddro— 4h f retder,(rddro+ ... (7a) 


tk L.—c, (0)—4k? j TAC (rad rod- KON (7 b) 
IU 
Similarlv, the expressions for the other coefficients become 


—2M,—2 [ Ses rod es — 4 f To ACa, (TdT t eee 


1 C) 


—3k M, =en,'(@)—4k? IN rode m.(redrot ... (Td) 


MALI =i euo) 2 [. [6i (rd—rAe(rd]drot « 
(re) 


2E L/—c'( »)— [^ Aei,(rddra—2k* |” [27 Ae: (rg) 
— To ACi (rd|drot... (7H 


—PMI- cuo) +2? f " [hes Cr — robes, (rg)]d rot. - 
2 : (72) 
-EMi =en I) [Aem rides 2k | [2roden,'(rd 


—réAcs(rg]dre-- ... (Th) 


It is interesting to examine these equations briefly with 
regard to the problem of one-degree-of-freedom oscillatory 
instability. «As was shown in reference 5, an oscillatory 
stability boundary is given by the equation .}f,=0. Hence, 
a wing pivoted about its leading edge can be neutrally stable 
if M;'—0. Such a condition arises only if the frequency is 
very low * for which case the above expansion for kJ’ applies. 
Hence, to a first order the stability boundary is given by the 
condition 


ee Cn) | Aen," rabd ro 


This leads to the interpretation of the indicial curves shown 
in figure 3; namely, an airfoil pivoted at its leading edge can 





Chord lengths traveled, T , 


fhea larra 


"ues 7 aa 


Le 


Oscillatory instability possible if 
€. (00) < fA tat, (To) dT, 


Ce a 


Chord lengths traveled, Te 


FIGURE 3.—Interpretation of singie-legree-of-freedoni stability boundary fn terms of 
indicial response curves. 


- A pulse function occurs only for the case }fo=-0. Its treatment in expressions such ss equations (5) is diséussed in a subsequent section. 
i That is, the moment of inertia Is large. Experimental verification of the existence of this type of instability can be found [n reference 10. 
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have an oscillatory instability if the steady-state value of 
Cm.” is less than the area between the steady-state value of 
Cm, and its indicial value when expressed in chord lengths 
traveled. 


PART H—METHODS FOR SOLVING UNSTEADY LIFT 
PROBLEMS 


BASIC EQUATIONS AND SOLUTIONS 


In terms of the perturbation velocity potential ¢, the wave 
equation can be written 


1 
—3 Pare (8) 


srt Sy t Vua — 
Gg 


where à, is the speed of pressure propagation, ¢’ is the time, 
and x,y,z are spatial coordinates. This equation applies to 
& flow field which is stationary at large distances from the 
disturbance region; furthermore, the coordinate system is 
stationary relative to the fluid infinitely distant so that, if 
& moving wing is being analyzed, the wing moves with 
respect to the z,y,2 axes. The last equation can be put in 
a more convenient form by introducing the notation 


t= ot" + 


so that the dimension of ¢ is length just as are the dimensions 
of the geometric variables z,y, and 2. Equation (8), 
together with this transformation, yields the canonical form 
of the wave equation 


x ey ss — Or: (9) 


and it is this form which will be considered. 

The first task 1s to study the relation between the motion 
of the wing and the coordinate system. As has already been 
mentioned, equation (9) is valid for a flow field produced by 
& wing moving relative to a coordinate system fixed with 
respect to the fluid at infinity. It is pertinent to consider 
the possibility of finding a transformation which will fix 
the origin of the coordinate system on the wing and, at the 
same time, retain the wave equation as the governing equa- 
tion of the flow. Certainly the first of these i is 
simple to fulfill if the second is neglected. 

The following transformation (known in Kisses theory 
as & Lorentz transformation) 


= a— Met 
J41— Mj 
n=y 
(10) 
(—2 
i— My 


where Mee, will satisfy both the above conditions. 


For Badio sippi that the wing is moving along the z 
axis with velocity V, Application of equations (10) to 
equation (9) makes the origin of the new system of axes also 
travel along the z axis with a velocity V, This is seen to 
be so since & is always zero when z= F =M. Hence, 
the £ axis is “fixed” on the wing. As to the second condition, 
a straightforward exercise in partial differentiation yields 


Pret Po F Pre = Pr (11) 


so that, in going from x,y,z to £,4,f space, the wave equation 
remains invariant; consequently, both the requirements 
mentioned have been fulfilled. 

It is instructive to consider briefly the consequences of 
applying the Lorentz transformation. Although the wave 
equation remains invariant, such physical quantities as 
length and pressure do not. For example, a wing with a 
chord c in the z,y,z space has, according to equations (10), 
a chord c//1— M in the &,4,t space. Furthermore, the 
loading coefficient which, on the basis of linearized theory, 
Is given in the z,y,z,! space by 


Ap 4 Og 


do VoM, Òt ud) 
becomes for the &,y,¢,7 space 
lE LL M (13) 
Yo V,J1— My AL, Or Ot 


If the wing motion is steady and there are no transient effects, 
equations (11) and (13) are independent of time and, to- 
gether with the resulting length transformations, become 


Pret Pay rr 


(14) 


These are immediately recognized to be Laplaee's equation 
for incompressible flow and the familiar Prandtl-Glauert 
compressibility corrections.® 

The preceding discussion has an important qualification, 
however, in the fact that the velocity of the moving co- 
ordinate system cannot exceed the speed of sound. A glanee 
at equations (10) serves to verify this statement since those 
equations show AM, must be less than 1 in order that and r 
be real for real z and é. In fact, it has been shown that there 
is no transformation which will fix the moving system of axes 
in a wing traveling at a uniform supersonic speed away from 
the original fixed axes and still keep the wave equation in- 
variant. _ Therefore, for analyzing a wing 1D. supersonic 


§ The equation Se a Me would first read £e (z—Afuo)/ J 1— Mọ where fo is a constant representing the ttime required for the motion to reach Ita steady state. However, ther coordi- 
nate can always be translated to any fixed position without affecting any of the equations for potential, loading, eto. Such a translation is assumed to have been made fn equations (14). 
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flight, it is necessary to abandon one of the two proposed 
requirements: Either the coordinate system cannot be fixed 
in the wing, or the field equation must be modified. The 
latter of these two alternatives has been studied by several 
authors (see references 2, 3, 4, 5) but it is the former which 
wil be considered in the present analysis. Further, since 
the axes cannot be made to travel as fast as the wing, they 
will not be made to move at all and equation (9) will be 
adopted throughout as the basic equation. 

Having decided upon the form of the partial differential 
equation, we must next establish the boundary conditions. 
For any given time these conditions sre similar to those 
studied in steady-state thin-airfoil problems; namely, either 
that the given slope of the wing surface is proportional to the 
vertical induced velocity ø, over the region occupied by the 
wing in the z=0 plane, or that the prescribed surface 
pressure is proportional to the timewise gradient o, in 
velocity potential over the same region. The addition of 
time simply means that this region moves about in the 
z=0 plane in conformity with the direction and velocity of 
the wing. 

The solution to equation (9), subject to the boundary 
conditions just mentioned, can be expressed by & formula 
which may be regarded as requiring either the evaluation of 
a double integral or the solution of a double-integral equation, 
depending upon whether & boundary-value problem of first 


or second kind is considered. This solution is known as . 


Kirchhoff’s formula. (See reference 1.) It may be written 
in a form convenient for aerodynamic applications as follows: 


1 1 
im -xJ m (z 2 


Ò , 
Oz ei; YiZis i—ry— 
Zt 


s 2 ru qa. 0 t) dS (15) 


where the A indicates the jump ( value on the upper surface 


minus value on the lower surface when applied to ¢ or se) 


of the function in passing through the z,=0 plane, 
r—4(r—zi-(y—y)-4G-—2), r= {rr t(y— yi? +23, 
and where the area of integration S, will be discussed in 
more detail later. 

The terms in the integrand of equation (15) can be 
shortened by introducing the following notation: 


1.0 _1 [òg 
Ta I: eii yv zut —n)-- 2 3z 2j (16a) 
AS. = Ò ox), Ys i— Da (2; Ór 2 eee 0, f— eet) 


MONTE | (16b) 


In this notation, the subseript r in equation (16a) means 
that r is to be held constant in the differentiation, and the 
prefix A obviates the necessity of indicating that the func- 


390 


tions considered are to be evaluated for z;—0, since it indi- 
cates that the difference of the values of the function across | 
the z—0 plane is to be taken. The right-hand side of 
equation (16a) can be recognized as a term representing & 
source located in the z, —0 plane, and the right-hand side of 
equation (16b) is seen to represent a doublet located in and 
with axis normal to the z,=0 plane. The brackets [] about 
the functions in equations (16a) and (16b) have & special 
meaning which is defined in the following way: If f is & 
function the value of which at a fixed point P depends upon 
the coordinates zi, yi, z;, § of a moving point Q, so that 


f-—f(zuysz-) 
[f] f Gy zt 7) 


then, 
- (17) 


where r is the distance from P to Q. As an example, con- 
sider the potential at a point P due to a moving source, 
the location of which at any time is Q. Then eg satisfies the 
condition just mentioned that it depends on the coordinates 
I; Yı 21, tof Q. The brackets []} indicate that the potential 
[e] depends not upon the source strength now at “time” t, 
but rather upon the source strength that existed “time” r 
ago. For convenience, [e] is referred to as the retarded 
value of ¢. 

The expression for a doublet (equation (16b)) is usually 
expanded as follows: 


E Ò = Ore Ó _ Orel 9e 
AE | x elaz .(: ro] O2 rs As, 
10 Ae? E) 
CE Fol , 0? 
du [Ag] Z +) orl z| 
Finally, equation (15) becomes 
ZA : L órg[ OAg 
eld i ind oe | c&elaz( Fo sola jes 
(19) 


The application of equation (19) awaits only a discussion of 
the area S, over which the integration is to be made. This 
discussion is important enough, however, to merit considera- 
tion in some detail and will be given in the following section. 


THE ACOUSTIC PLAN FORM 


Suppose that a line of sources is placed along the y, axis 
and that the strength of these sources is zero for t<0. At 
¿=0 they are “turned on" and, at the same time, start 
moving along the negative x, axis with the velocity Y. 
After time 4# has passed, the source line has traveled & 
distance M as shown in figure 4 (which is drawn for the 
case Jf,>1). Suppose next that there are two sensing 
elements, or detectors, placed at the point P(x,y) located 
somewhere ahead of the y; axis; one of these detectors is 
responsive to light and the other to sound. Now, the light 


$ The z=0 plane is assumed to be the “plane of the wing’’; that is, If the angle of attack were zero and the wing had no thickness It would lle entirely in the z=0 plane. 


T Quotes are used around the word time since the dimension oft Is actually length, not time. It [s convenient, however, to refer ta fas "time," and, since the actual valueoftime is simply f 


t divided by the constant ae thts shouki cause no confusion. 
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FIGURE 4.—Acoustic wave pattern for a moving line of sources. 


detector at any given time will show the sources lying in a 
straight line just as they would appear visually at each 
particular instant. The situation is entirely different, 
however, for the sound detector. First it is necessary to 
understand the nature of a spherical sound wave. Such a 
wave travels outward from its origin at a velocity @, so that 
in the time ¢’ it has traveled a distance t. Before the wave 
reaches a point, the point is completely unaware (unaware 
is used in the sense that an instrument will record no change 
in any of the physical properties of the air at the point in 
question) of its existence and, further, after the wave has 
passed, the point remains subsequently unaware of its 
existence. Hence the only points disturbed by the wave 
are those momentarily on the spherical surface itself. (In 
this connection see reference 1, pp. 1-3.) 

The sound detector, therefore, can only “hear” sources 
which are so located that their spherical sound waves are 
just, at the given instant of time, reaching the detector. 
The locus of all the points which, at a time ! ago, emitted 
sound waves that are just now reaching the point P(z,y) is 
itself a sphere and for convenience this sphere will be referred 
to as an “inverse sound wave."? The traces of these inverse 
sound waves in the z=0 plane are drawn in figure 4 as 
concentric circles about the point P(z,y). The intersection 
of an inverse sound wave of radius £—- with the line repre- 
senting the position of the sources at a time !'—7' ago gives 
the position of the sources which are just now signaling their 
presence to the sound detector at P. For example, when the 


source line started, it was lying along the y, axis. With 
reference to the present time ¢ this was aj removed. Hence 
the intersection of the circle about P of radius £ with the y, 
axis fixes the two points A and A’, the spherical sound waves 
of which are now reaching P. A continuation of this process 
yields, for the locus of all points from which waves emanated 
that are just now touching the point P, the part of the 
ellipse. shown in figure 4. 

In the sense that the light detector, because of the very 
large velocity of light, is "seeing" a straight line of sources, 
the sound detector, because of the relatively slow velocity 
of sound, is “hearing” an elliptic line of sources. Extending 
this concept to include a sheet of sources distributed over 
the surface of the wing, one can refer to the outline of that 
part of the wing which generates disturbances which can be 
measured by the light detector as the plan form (i. e., the 
visual plan form), and to the outline of that part of the wing 
which affects the sound detector as the acoustie plan form. 
In a mathematical sense, the acoustic plan form is the area 
S, over which the integration of equation (19) is to be made. 

The equation for the acoustic plan form can be formulated 
by means of the two equations 


(z—z + (y—y)! + 2— (t— 7) (20) 
Sy, 2, 7) —0 (21) 


where z, y, z are the coordinates of the point at which tho 


induced effects are to be measured; z;, y;, the variable points 
of the sources; t, "time" now; and í(—7, “time” ago. Equa- 
tion (20) is that of the inverse sound waves and equation 
(21) represents the position of the visual plan form at a 
"time" z. It is necessary to include the region behind the 
wing covered by the vortex wake as part of the visual plan 
form. In case the vorticity in the wake vanishes, as in the 
thickness problem, the wake may still be considered as part 
of the visual plan form, but the strength of the source-doublet 
distribution over that part of the acoustic plan form corre- 
sponding to the wake will vanish. If @—2?<00, the acoustic 
plan form does not exist since, for such a case, no source 
on the zı, y; plane has had time to transmit its effect to the 
point z, y, z. On the other hand, if the circle (in the 2, yi 
plane) given by the equation 


(x—2zn)'--(gy—y)^-— 


les entirely within the area occupied by the visual plan 
form of the wing at the beginning of the motion, the acoustic 
plan form is just this circle itself. For any other situation 
the acoustic plan form is formed in part, or in whole, by the 
curve found by eliminating r from equations (20) and (21) 
and in part, or not at all, by an are of the circle (z—2;)*-- 

(y—13,)-8—2, One of the principal advantages of the 
acoustic-plan-form concept arises in the study of problems 
involving source or doublet distributions having constant 
| strength. In such cases, the retarded values of the potential 


= (?— z2 


‘Tha inverse sound weve has for its analogue in steady lifting-surface theory the Mach forecone. In that theory a disturbance outside the Mach furecone cannot affect the values of 
the induced velocities at the point where they are being measured, Similarly, in the present study, a source located outside the inverse sound wave of radius č cannot affect the valucs of 


any measurement made at the pobit P. 
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and its gradient appearing in equation (19) are constant 
and can be taken outside the integral signs. The problem 
is thereby reduced to the integration of a simple geometric 
variable over S,. . 

A few examples will serve to fix the idea of the acoustic 
plan form. Consider first a two-dimensional, unswept wing 
moving at a constant supersonic speed in the negative rı 
direction. At time zero the leading edge of the wing was 
along the y; axis and now, at time equal to //, the wing has 
moved so that the leading edge coincides with the line 
D, —JM,. Choose three points that are now lying on the 
wing. Let one point have its z coordinate in the range 
e—4U—232z24U—z* (where e is the chord of the wing), 
the second in the interval f@—z*>r>—,f—2?, and the 
ihird in the range —£—2?»z2—AM. Designating these 
points by P,,P,, and P; (see fig. 5), it can be shown that their 
acoustic plan forms are, respectively, a complete circle, 
a part circle and part ellipse, and a complete ellipse. The 
points P are at the centers of the circles and at focal points 
of the ellipses. Since, moreover, the circular plan form about 
P, receives no signals from sources on the leading or trailing 
edge, conditions at P, are consequently completely inde- 
pendent of the actual (visual) plan form of the wing. The 
elliptical plan form about P;, on the other hand, depends 
entirely on the shape of the leading edge; and finally the 
mixed plan form about P, is in certain regions (the circular 
portion) independent of the leading edge, and in other 
regions (the elliptic portion) entirely dependent upon it. 
Since the wing is traveling at supersonic speeds, the trailing 
edge and vortex wake can have no effect on the measure- 
ments taken on the wing and, in the same way, a point ahead 
of the wing leading edge, P, in figure 5, is undisturbed. 

Next consider a wing moving at a constant subsonic speed 
in the negative x, direction. As before, the leading edge 
was on the y; axis at é’=0 and has traveled a distance — Jfa. 
Choose now three points P,, Pa, and P; on the wimg and 
imaffected by the wing tips. The acoustic plan forms for 
these points are combinations of circles and hyperbolas as 
contrasted with the circle-ellipse combination m the super- 
sonic case. Just as in the supersonic case, however, there 








zy 


———— Wing at fime zero 
Wing now 


FiGURE 5.—Acoustie plan forms for points on wing traveling at supersonic speed. 
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is a certain region represented by P; in which the acoustic 
plan form is a complete circle and is independent of the 
visual shape of the wing (see fig. 6). Point P, is surrounded 
by a plan form which is part hyperbolic and part circular, 
the point itself being the center of the circle and the focus 
of the hyperbola. Point P; is a limiting value of Pa; 


it lies on the leading edge of the wing and the hyperbolic 


sides of its plan form have degenerated into straight lines. 
Finally, P, lies ahead of the wing; its plan form is still a 
combination of a hyperbola and a circle, but P, is now the 
focal point lying ahead of the hyperbolie branch used. 
Figure 6 was constructed so that the portion of the visual 
plan form behind the trailing edge had no effect on the 
potential at the various points Pi, ete. If these points had 
been chosen at positions where the wake could signal its 
effect, one of two acoustic configurations would result. 
First, if the wing is symmetric about the 2—0 plane, no lift 
is developed and the vortieity in the wake is zero so that the 
visual plan form need not include the wake, but effectively 
ends at the trailing edge. In this case, the leading edge of 
the acoustic plan form is then determined as before, while 


its modified trailing edge may be made up, in part, of circular _ 


ares formed by the primary wave and, in part, by an arc of 
the hyperbola formed by the (acoustic) intersection of the 
straight visual trailing edge with the primary wave (such an 
arc being identical with the leading edge of the acoustic plan 
form but displaced backwards). On the other hand, if the 
wing has no thickness but is inclined to the free stream, it 
develops lift and the vortieity in the wake does not vanish; 
the acoustic plan form has a trailing edge made up entirely 
of an are of the primary inverse sound wave. The space 
between this arc and the acoustic trace of the visual trailing 
edge is covered by a sheet of doublets, the strength of which 


is determined by the vorticity distribution of the vortex | 


wake. 

It is interesting to notice the conversion of terminology 
which arises in the analysis of unsteady lift problems. In 
the study of steady-state wings, it is customary (because of 
the nature of the governing partial differential equation) to 
speak of the subsonic problems as elliptic and the supersonic 


E 


x 


————— Wing af time zero 
Wing now 


Figure 6.—Acoustic plan forms for points on wing traveling at subsonic speed. 
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problems as hyperbolic. Yet the acoustic plan forms just 
presented involved ellipses for the supersonic wing and 
hyperbolas for the subsonic case. 

To complete the remark, it can be observed that when the 
velocity of the wing is sonic the steady-state partial differ- 
ential equation becomes parabolic and, in this case, the 
acoustic plan form of a straight-edged wing also involves 
parabolas. When the leading edge is linear and normal to 
the stream direction, the eccentricity of the conie sections 
hounding the acoustic plan form for a point on the wing is 
equal to 1/Af. From this relation it is apparent that for 
A, less than, equal to, and greater than 1 the sections are, 
respectively, hyperbolas, parabolas, and ellipses. As might 
be presumed, the value of the eccentricity satisfies simple 
sweep theory so that, for an infinitely long straight leading 
edge, the eccentricity of the acoustic plan form is 1/44 cos A 
where A is the angle of sweepback. The principal axes of 
the conie sections are always normal and parallel to straight 
leading edges. 


HOMOGENEOUS BOUNDARY-VALUE PROBLEMS 


Kirchhoff’s solution to the wave equation can be applied 
to arbitrary wing plan forms undergoing arbitrary maneuvers. 
The boundary values for such general problems, however, 
usually lead to the development of double integral equations 
which are difficult to solve. As is usual in such cases, there 
are many special types of plan forms and maneuvers which 
lead to boundary-value problems that are simpler to analyze. 
An important class of these simplified problems is that 
arising from homogeneous boundary conditions. 


Let e(z,y,z,7 be a solution to equation (9). In certain 


special cases this can be written p= (t)"q, (? 5 i) in which 


case e is called a homogeneous function of degree n. The 
number of variables affecting g is only three as compared 
to the four which are necessary to determine e. If, there- 
fore, a partial differential equation can be set up for e, it 
wil contain one less mathematical "dimension" than the 
equation for e. Following this observation it is necessary 
to proceed in two directions; one to find the partial differential 
equation for g, and the other to find the physical problem 
and consequent boundary values leading to a homogeneous 
flow field. The latter path will be first explored. 

First, consider an example of & homogeneous boundary- 
value problem. Suppose that & rectangular flat plate starts 
suddenly from rest and moves forward at an angle of attack 
at a supersonic Mach number M, At “time” f, the initial 
spherical wave generated by the forward right-hand corner 
has traveled outward to & radius ¢, and, at “time” 25, to a 
radius 26. Figure 7 indicates the traces of these spheres in 
the z=0 plane together with the original and present position 
of the wing leading edge. Let the points P, and P, be 
located on the same rays through the origin of the circles 
and the wing corners. "The problem is to find the pressures 
at P, and P3. 

It is apparent that, 1f every dimension in the figure involv- 
ing P, is divided by 21, and every dimension in the figure 
involving P, is divided by &, the two figures will be similar 
in every respect and point P, will coincide with point P}. 
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Friavgz2 7.—Geometric relationship for homogeneous flow. 


Since the vertical velocity w is constant over the plan form, 
& simple change in scale has made the boundary conditions 
for both problems identical. But this means that the solu- 
tions at P, and P; are identical since the wave equation is 
invariant to change in scale. Hence, in regions of a rec- 
tangular wing unaffected by the waves from ihe trailing 
edge, the pressure can be written 

Ap Ap (5 yg 

p vw a 
and the pressure is a homogeneous function of degree zero. 


A generalization of this example is contained in the following 
statement: 

(1) The pressure in any region affected by only two 
intersecting edges of a straight-sided flat plate traveling 
at a uniform subsonic or supersonic speed is homogene- 
ous and of degree zero (1. e., satisfies equation (22)). 


TWO- AND THREE-DIMENSIONAL UNSTEADY LIFT PROBLEMS IN HIGH-SPEED FLIGHT 


Consider as another example the case of a flat rectangular 
wing traveling forward at a subsonic or supersonic speed and 
rolling about one edge, taken to be coincident with the z axis. 
The argument follows the same lines as before, and again a 
change in seale, proportional to the time, makes the geometry 
of the wing-wave combinations identical (in regions affected 
by only two intersecting edges} for different times. The 
boundary values over the wings will not be the same, how- 
ever, unless the slope of the w, distribution is adjusted in 
each case. But wy can be adjusted by reducing it an amount 
proportional to the distance from the axis of rotation. The 
boundary-value problems are then similar for different values 
of time. Finally, therefore, the pressure can be written 


AD isis esp = SL) 
Qo (£, Y, 2,0) — Qo (> tt 


which is a homogeneous function of degree one. 
ization of this example is expressed as follows: 


(23) 
A general- 


(2) The pressure in any region affected by only two 
intersecting edges of a straight-sided flat plate traveling 
at a uniform subsonic or supersonic speed and rotating 
at a constant rate of pitch or roll is homogeneous and of 
degree one (1. e., satisfies equation (23)). 


It should be noted that both (1) and (2) are equally true 
for the steady-state case when all transient effects have 
disappeared. In supersonic wing theory they lead to conical 
and quasi-conical flows, respectively while in the subsonic 
case they lead to flows about wings having infinite chordwise 
extent. In general, homogeneous flow occurs when the 
boundary conditions after & change in scale are proportional 
to their original values. 

Consider next the modification of the basic partial differ- 
ential equation (equation (9)) under the assumption that the 
flow is homogeneous. If the pressure is given by a function 
that is homogeneous and of degree zero, then, by equation 
(12), the velocity-potential function will be homogeneous 


and of degree one. If the notation 
Y ? Z 
E Io t= yo, p“ 
(24) 
e, y;?, t}=tP(2Xa, Yo: Zo) 
Is used, then equation (9) becomes 
(1— rð Prazo t (1— yo) Py po V (1— 20") Pr uo — 22040 Pry — 
2 Lozo Prg 2Jo2o Prg (25) 


and a linear partial differential equation with three inde- 
pendent variables is therefore obtained. 

In the general theory of partial differentia] equations of 
second order, the character of an equation is determined from 
the geometric nature of a related quadric surface. The char- 
acter of equation (25) can be shown from such considerations 
to depend on the sign of the expression 1—x*—y’— zy. It 
is immediately apparent, however, that within the unit 
sphere in the 29,¥0,20 space the sign of 1—2z—yj—z4 is 
everywhere positive and outside the sign is everywhere 
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negative. It follows that outside the unit sphere equation 
(25) is hyperbolic and inside the unit sphere it is elliptic. 

The character of equation (25) is of particular interest 
since the difficulties inherent in the determination of the 
solutions can be estimated without actually obtaining the 
solutions. For example, consider. the two configurations 
shown in figure 8. These wings started moving at £—0 with 
the foremost portion of their leading edges on the ys axis and 
have by now traveled forward at a supersonic speed to attain 
the positions represented by the figure, the unit circle being 
in each case the trace of the primary wave from the vertex on 
the z—0 plane. Outside the unit sphere, the governing 
equation is hyperbolie and the behavior of the flow is similar 
to that in steady-state supersonic-wing problems. Inside 
the unit sphere, on the other hand, the character of equation 
(25) is elliptic. 

It is instructive to notice that this entire development has 
a direct analogue in the study of three-dimensional, steady- 
state, supersonic wings. In that case the original equation 
is the three-dimensional wave equation 


Qz— y Pu = 9 (26) 
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FIGURE $.—Homogeneous flow flelds. 
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By considering the velocity potential to be homogeneous and 
of degree one, Busemann in reference 11 was able to introduce 
the transformations 


yw 
g(z,9,2)— rb (yo, 20) 
and transform equation (26) to the form 


(1— yo) Brog (1 — 297) $4, — 2yo2 45, ,,— 0 (27) 


which is the two-dimensional form of equation (25). Flows 
governed by equation (27) have become known as conical 
flows. A study of equation (27) shows it to be elliptic inside 
and hyperbolic outside the unit cirele. In this case, however, 
the equation has only two independent variables so .that 
once the equation has been transformed (by means of the 
Tschapligin transformation) to the two-dimensional form of 
Laplace’s equation solutions are not difficult to find. 

The simplification of the four-dimensional theory brought 
out by the introduction of homogeneous flow was more ap- 
parent than real since the resulting partial differential equa- 
tjon, although containing one less dimension, was unwieldy. 


BOUNDARY-VALUE PROBLEMS INVOLVING NONINTERACTING SURFACES 


Another class of wing problems which is simplified both 
in theory and in practice by reasoning from physical knowl- 
edge of the flow behavior is that in which the wing has a 
supersonic edge (i. e, an edge which. is traveling with a 
supersonic normal component of velocity). 

When the acoustic plan form is affected only by & super- 
sonic edge, it is apparent that the flow on the upper surface 
of the wing is independent of that on the lower surface. 
Hence the solution to such problems can always be written 
in terms only of sources as follows: 


EE 1[0e 


where 99/02; = wu(21,91) = Vou (21,51), ^y being the local slope 
of the surface in the direction of Vo. Since the equation (28) 
is equally valid for symmetrical nonlifting surfaces and lifting 
plates, its value and simplicity are evident. 

If the wing plan form is further specialized by having not 
only supersonic leading edges, but also having a straight 
trailing edge perpendicular to the direction of motion, addi- 
tional simplifications can be used.® Consider, for example, 
the two-dimensional wing (a) in figure 9. Let this wing have 
an angle of attack a(é) which varies with time im an arbitrary 
manner. There results from such an angle-of-attack varia- 
tion & certain lift which also varies with time. Hence, if 
L* represents the total lift on an airfoil of very high aspect 
ratio and c;,* represents the section lift-curve slope, then !? 


ere CE CTT 
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Next it is clear by reason of symmetry that the total lift on 
wings (b) and (c) in the figure are equal. Then, since the 


just half of that on wing (a). 
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FIGURE 9.—Dovelopmont of elemental lifting strip. 


analysis is based on a linear partial differential equation, by 
superposition principles the total lift on wing (b) or (c) is 
In another sense, the lift 
coefficient for the whole wing based on the deflected area is 
the same for all three cases. A suitable superposition of 
wings (a), (b), and (c) will give wing (d), which then has 
the same lift coefficient based on deflected area. Finally, 
because of the supersonic stream, wing (e) can be obtained 
from (d), hence it also has the lift coefficient common to the 
other wings. It is, of course, necessary that the variation 
of æ with the time be the same in each case. 

The preceding process can be extended one step farther 
to the development of the lift due to a single deflected 
element. By considering figure 10, it ean be seen that. 


L* ce;.*(e,t)—(e—Ac)e; *(c—Ae, è) 
ert Mori AD caesa Sod LEE V. 
Joa Ac 


where AS is the area of the deflected element and c is the 
distance_from the centroid of AS to the trailing edge. By 
the usual limiting process the latter equation becomes 

L* 


=< [e e;,* (e, 0] dS 


t The following method simply extends, to Include the effects of unsteady motion, a theorem given by Lagerstrom and Van Dyke. (Bce reference 12.) 
B The asterisks on quantities indicate that two-dimensional] values are taken, or that a high-espect-ratio wing is considered and tip effects are neglected. 
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Fr;vEE 10.—Development of lifting element. 


Finally, if a wing is composed of a distribution of these 
elements, then, for a coordinate system centered at the 
apex of the wing leading edge, there results 


I 
Cag h [ar Meze coz idydr (29) 
where c 1s the maximum chord (see fig. 11). In the develop- 


ment of equation (29) each element is assumed to have the 
same variation of motion with time. 
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Firer 11. — Example of use of lifting element. 


Notice that when the wing is a flat plate flying at a steady 
speed so that all transient effects have disappeared, c; * is 


independent of ¢ and of the chord length, being, in fact, 
equal to 4/8. Then equation (29) becomes 


4 4a 
Cs [, f «8 zv 


where @ is the average angle of attack of the wing. This 
result has already been obtained in reference 12. When 
aír, y) is independent of x (as for a flat wing sinking or 
rolling), equation (29) becomes 


VE : . 
=g f CO € €; (c, D)dy (80) 


where c is the local chord which is, in general, a function of y. 
Equation (30) simply indicates that longitudinal strip theory 
is exact for calculating the total lift on such wings. 

Finally, notice that the calculation of the unsteady lift on 
three-dimensional wings with supersonic leading edges and 
straight trailing edges perpendicular to the free-stream 
direction has been reduced to an integration involving the 
relatively simple results for & two-dimensional wing under- 
guing the same unsteady motion. For example, the lift 
on a flat unyawed triangular wing with supersonic leading 


edges rising and sinking with & harmonic motion can be 
computed from & single integration of the results presented 
in reference 5. 

There is another simplified method for obtaining the total 
lift and moment on a wing with all supersonic edges and a 
straight trailing edge. In this presentation it will be assumed 
that the trailing edge is normal to the free stream. How- 
ever, since the wave equation is invariant to a rotation, it 
will be apparent that the solution can be generalized to in- 
clude a straight, supersonic trailing edge yawed with respect 
to the free stream. 

Consider equation (9) and integrate each term with respect 
to y between the limits minus and plus infinity." There 
results the equation 


sd (^ gayt [^ kay- (7 5 dy—o 


If y yi(z, 2, 0 and y=y,(z,2,f) are the equations of the Mach 
waves streaming back from the leading edges on the left and 
right sides of the wing, respectively (see fig. 12), then, since 
e is continuous across these waves but gz ey, v, and e, 
&re not, 
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where u, and v; are the values of 4 on the interior faces of 
the right and left Mach waves, respectively, and 
Fr oe 
n dy? ,dy—9,— 


Values of the terms involving e, and e, are similar to those 
involving e; so that finally, if 
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FiGCRE I2.—Forward portion of Mach wave system for a supersontc-edged, triangukir wing. 


\t The basic idea for this solution was given by P. A. Lagerstrom in his lectures at the California Institute of Technology. 
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The terms enclosed within the parentheses in the last 
equation combine internally so that each is zero. For the 
case of interest here, this is not difficult to show. Consider, 
for example, the Mach wave streaming from the right edge. 
Then, since the equation of this wave is 


y,—-—z4g^m:—1 +mz +m Mot 


and the value of the potential is the steady-state two- 
dimensional value given by the expression 


X m(rd-Md)—y—ed8*m—1 
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and this is identically zero.? Finally, therefore, equation (9) 
has been reduced in terms of equation (31) to 


bi: — 5, — 0$, —0 (33) 


The boundary condition in terms of ® for a triangular 
wing with supersonic edges is given by 


Op EE: Vy ER "Oe 
a yi edy- f; Oz 


where w, is the vertical induced velocity in the plane of the 
triangular wing. The derivative with respect to z can be 
carried through the integral sign because the extra terms in- 
volving the value of » at y; and y, vanish. In fact, since 
equation (34) applies to the 2—0 plane, the limits y; and y. 
can be replaced with the expressions for the left and right 
leading edges of the triangular wing, respectively. The 
boundary condition expressed by equation (34), used in con- 
junction with equation (33), suggests a problem exactly like 
those posed by lifting surfaces in steady-state wing theory; 
in fact, the problem of & wing tip of specified camber in a free 
stream at Mach number4/2. Figure 13 shows a lifting sur- 
face in the z,t plane. The solution for the potential & in the 
steady-state problem can be written 


Ob 
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dy= | w dy (34) 
raf Jt 





where ¢ is the area on the wing plan form in figure 13 that is 
included in the forecone (t—4 Y = (z—2j)*. 
Now from equation (12) we have 


AP quus. 92 
LL go dy— VAM, t s=6 
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FIGURE 13.—Region of rt plane in which boundary conditions for P aro known. 


where s is the local semispan of the triangular wing. This 
equation shows how the solution to the lifting surface prob- 
lem in æ will aid in the solution of the unsteady problem, for 
the unsteady lift on the triangular wing is given by 


u Í ĉo & Ap e 4 €g oo 
C. Í, dz [^ d dy—gatv. |, ot 


where S is the area of the triangular wing. It is therefore 


dz 





seen to be convenient to evaluate the quantity 2 , given 
3m0 
by 





Ob 
OF (35) 
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The pitching-moment coefficient for the triangular wing in 
unsteady motion can be evaluated similarly. Specific appli- 
cations òf this method will be found in & subsequent section. 





£-0 


TWO-DIMENSIONAL BOUNDARY-Y ALUE PROBLEMS 


The simplification brought about when the flow is inde- 
pendent of one dimension is obvious. In such cases, the 
three-dimensional wave equation (8) reduces immediately to 
the two-dimensional wave equation. Typical examples of 
this type of problem can be constructed by considering flat 
plates which start suddenly at £—0 and travel thereafter at 
constant supersonic velocities. Two examples, one a corner 
of a rectangular wing and the other a triangular wing, are 
shown in figure 14. After time £—0, the edges of the wings 
send out cylindrical waves and the outer boundaries of these 
waves at time t are shown as dashed lines parallel to the edges 
in question, Since points in regions 1 and 2 are affected only 
by a single edge, the wave phenomena in these regions are 
cylindrical, and the physical quantities are in both cases 
independent of distance parallel to the edge which acts as 
their generator. Hence, the flow field in these regions is two- 
dimensional. (Region 3, incidentally, is independent of 
distance in both z and y directions and is, therefore, one- 
dimensional.) 


13 It is not necessary to perform a direct calculation m order to prove the above result for arbitrary plan forms. The terms m parentheses in equation (82) represent the directional deriva- 
tive of the velocity potential teken along the so-called "oonormal'' of the foremost disturbance surface. Since e is constant on the surface, and since the conormal lies along the surface, the 


terms in parentheses are zero. 
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Ficat 14.—Cases of two-dimensional flow flelds. 


Solutions to the two-dimensional unsteady problems are 
sometimes especially easy to find because of the analogy they 
have with three-dimensional, steady-state, lifting-surface 
problems. For example, consider an infinitely long unyawed 
wing which starts from rest and travels forward at a velocity 
Ve which may or may not be a function of time. The trace 
of this wing in the y, t plane is like that shown in figure 15. 
(In the figure shown, the wing velocity is varying and is 
always less than the speed of sound.) The boundary condi- 
tions are that ¢, is specified over the shaded area and the 
loading Ag, is zero everywhere except within the shaded 
area. But if z is replaced by y and t by x, these boundary 
conditions are exactly the same as those for a plate of known 
camber and angle of attack, with a plan form as indicated 
by the shaded area, placed in a free stream directed along 
the positive z axis at a Mach number equal to 42. The 
solution for the one problem may be used, therefore, as a 
solution to the other with onlv a change in notation. 
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FiG'vgE 1§.—Decelerating wing in x plane. 
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BOUNDARY CONDITIONS FOR VERY SLENDER WINGS 


When the wing plan form is slender in the sense that its 
length in the streamwise direction is large compared to its 
span, an estimation of the loading on it can be obtained by 
neglecting in the partial differential equation the gradient 
of the induced velocity component in the stream direction. 
Thus, if the wing is moving in the negative z direction, 
equation (9) reduces to 

£y t 9791 (36) 
which is again the wave equation but in two space dimen- 
sions. Since equation (36) is independent of r, study can be 
made independently in each plane z—constant. This is an 
extension of steady-state slender wing theory, see, e. g., 
reference 13. Figure 16 (a) shows a typical section in the yt 
plane. If the wing is a flat plate at & constant angle of attack, 





Section 


(b) 


{a) Axes fixed relative to still afr at Infinity. 
(b) Axes fixed on wing. 


FIGURE 16, —£lender wing in unsteady flow. 
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the value of e, over the part of the yt plane occupied by the 
wing is & constant and the jump in e across the vortex wake 
must be consistent with its value at the wing trailing edge. 
The analogy with three-dimensional, steady-state, super- 
sonic, lifting surface theory is apparent. 

There is another way by means of which the effects of 
unsteady motion on slender wing pressures can be estimated. 
If instead of using the stationary zyzt coordinate system, the 
reference axes are fixed on the wing by the simple set of 
transformations 


z,—izd-Mst 
t =t 
yı =y 
21772 


equation (9) becomes 


(1—M?) ex, —2Mo 9: + ey T qi = 01 
Again, if the induced velocity components in the stream 
direction are neglected, the simplified equation 


results. The latter equation is identical in form to equation 
(36). However, now the axes are fixed on the wing and a 
typical section in the ¥,f; plane is similar to that. shown in 
figure 16 (b). In this case, a flat plate wing is represented by 
& constant value of e, over the entire shaded area in the figure 


and elsewhere ¢ must be continuous.” 


TWO-DIMENSIONAL UNSTEADY INCOMPRESSIBLE FLOW 


The analogy between two-dimensional unsteady and three- 
dimensional steady flow includes the case of a two-dimen- 
sional, unsteady, incompressible flow field the analog of 
which is a three-dimensional, steady flow field having a 
free-stream Mach number equal to 1. This can be demon- 
strated by inspecting equation (8). Since the flow 1s two- 
dimensional and since for an incompressible medium the 
speed of sound a, is infinite, the basic equation governing the 
flow can be written 


Pr t P=) 


It must be remembered, however, that time still appears in 
the boundary conditions and in the equation for the loading 
coefficient which, according to equation (12), can be written 


Hence, the basic partial differential equation and the expres- 
sion for the loading coefficient are the same, except for a 
change in notation, as those governing three-dimensional, 
steady-state problems when 44,=1. 

If a two-dimensional wing in an incompressible fluid 
starts from rest and travels forward at a speed Vo, the trace 
of the wing is as shown in figure 17. The essential difference 
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FiGURE 17.—T wo-dimenslonal wing in rë plane; Mo» 0. 


between this problem and the more general case of two- 
dimensional compressible flow lies in the fact that in this 
case the traces of the characteristic cones are normal to the 
t’ axis. The boundary conditions are therefore satisfied 
along lateral strips and, in lifting-surface terminology, the 
analysis corresponds to slender-wing theory. These latter 
methods are well established and in reference 13, for example, 
the manner in which the Kutta condition is impósed is dis- 
cussed in some detail. The trailing vortex sheet for the 
lifting wing has the same distribution of vorticity that exists 
behind the unsteady two-dimensional airfoil and the rolling 
up of the vortex sheet can be studied from either standpoint. 


PART IIIL—- SOME APPLICATIONS OF THE METHOD 
STARTING LIFT OF A WING 


One. of the simplest and yet most general resplts which 
can be derived on the basis of the present theory is the initial 
value of pressure on a wing surface starting suddenly from 
rest with a velocity Fo. The discussion will be made for a 
wing without thickness although the method will be seen 
to apply to the thickness case as well. 

Consider a surface with a plan form as indicatéd in figure 
18. The acoustic plan form of a point P(x,y) on the surface 
is a small circle of radius £. Since no point on the wing out- 
side this circle can influence the pressure at P, the upper 
surface is independent of the lower surface, except for a 
band of width ¢ around the edge of the wing. It is, there- 
fore, evident that the boundary-value problem to be solved 
has been treated in the section Boundary-Value Problems 
Involving Noninteracting Surfaces. The solution follows 
directly from equation (28) and can be written 


“as Ja) sns. | 58 


u The analogous problem in steady-state wing theory is that of a low-aspect-ratto, rectangular, flat plate in a free stream having a Mach number equal to 4/2. 
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FiGCRE 18.—Acoustie plan form at start of motion. 


Using a polar coordinate system defined by 


£—4,-r cos 8 

y—3y,—r sin 0 

da, dyy—rdrd6 
there results * 


NETUS (^as "ar — 
= 2r Is d6 | dr— toV, y) 
so that 
(22 E: van — Aw, y) 
. Qo7 t-o VoM, t-0 VoM, 


If æ= —w,(z,y)/Vo is the local slope of the wing, the expres- 
sion for load coefficient becomes 


(= Aer. y) 
Qo /t«-0 Af, 


The starting value of lift coefficient can, therefore, be. 
written 
io 


CL— M. 


where æ is the average angle of attack of the surface defined by 


z d due 
=g) f dx dy 


S being the area of the wing plan form. 
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| INDICIAL FUNCTIONS FOR A TWO-DIMENSIONAL FLAT PLATE, M4-0 


The basic partial differential equation governing an 
unsteady flow feld in the two-dimensional incompressible 
case is Laplace's equation in two dimensions 


Pr F Zu =l | (37) 


where x is distance along the chord and z is height above the 
plane of the wing. It was also pointed out in Part I that — 
the boundary-value problems arising in the study of the 
unsteady, two-dimensional plate flying in an incompressible 
medium were directly analogous to those which are studied 


in three-dimensional, steady-state, lifting-surface theory - 


under the classification “‘slender-wing theory." This analogy 
is useful since well-established concepts in one field can be 
immediately carried over into the other. It should be 
mentioned, however, that the subsequent treatment of the 
incompressible case is not intended to be an improvement 
on Wagner's original derivation (see reference 14) but 
rather it is a rederivation along lines that will be used later 
in the analysis of the compressible case. 

The initial pulse.— The first analogy with slender-wing 
theory which will be used concerns the initial pulses that 
occur in the values of Iift and pitching moment. It is a 
well-known result (reference 15) that the total lift, as given 
by slender-wing theory, on the wing shown in figure 19 (a) 
is & function only of the maximum span and the value of 
wy, along the section of maximum span (section AA). It is, 
therefore, independent of the wing twist and leading-edge 
shape ahead of section AA. This concept has been extended 
in slender-wing theory to the extreme case shown in figure 
19 (b) of a rectangular wing. The lift on such a wing is 
concentrated entirely along the leading edge and is a function 
only of the span of that edge and the value of w, there. By 


: the analogy existing between the two theories, therefore, it 


is evident that the solution to the indicial problems in two- 
dimensional, incompressible, unsteady flow (fig. 19 (c)) will 
contain a pulse at ¢’=0. 

The evaluation of this pulse will be treated briefly. A 
solution to equation (37) for the vertical induced velocity 
in the z=0 plane can be written in terms of the jump in u 
across the z=0 plane (see reference 13); thus, for the shaded 
area in figure 20 this is 


1 f^ Au(rj 


2r f a £—1 


dz, (88) 





w(x) = 


The general inversion of equation (38) can be written 





A 
A = Ã— 
u(x) sd zEaymz 
L2 1. A (Li) Fa EGXÉL. x 1 
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where 
A- [' Audz 


Wa (E, 
4 The mean value theorem gives 9**— gn ff n where £ and q He somewhere in Se. Hence, as £ approaches 0, £ and y approach r and p, respectively. 
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FIGURE 19,—Lift pulses In slender-wIng theory. 


In the present case A is zero, since Ag is zero at z— —a and 
3 —b, and an integration of both sides of equation (39) with 
respect to z between the limits b and z gives 


vor =a) (ab) az — z) dz. 


V (z —21) (a+ b) 


Adoption of the notation 


x(t’, zj)-- constant = Í wt’, zi)dzi 
and integration by parts leads to the equation (since 


BENE e e 
s c NEED vanishes if 575 17» — a) 


dus 2J(b—z)(a-- ) —*t',x)dz 


(40) 


T -a(z —2)4(5 — xi) lat zi) 


The loading can now be determined by using equation (12) 
(and GUIBCer Hn uig with respect to f rather than f£) 

T 4 òg 

— Vg ot 

If the shaded area in d 20 is allowed to vanish, all the 

loading accumulates along the x axis in the region 0S zSc. 

Therefore, the integral of the loading with respect to ¢ over 

the shaded portion must be considered. The final result for 

the pulse loading (Ap/qo), at t’=0 can be expressed in terms 

of the à function (see list of symbols) as 


An 45(t’) i, aN E WO, z tide, —— 
a avg ea ?(z—2)4—2)no 


The boundary : conditions for the sinking and pitehing wing. 
given by equations (2) and (3), when inserted into equation 
A), 


(41), viele 
A 
T V ve us 
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re 


(42) 
zr) 


After integration, the pulse values for lift and pitehing 


moment may also be obtained. Hence, 
(c1) y at) 
(Cna J=- 3T &(t^) 
MEM (43) 
(Ct, =T Xt) 
F 9 £ 
(m, h= 64 T. t ) |. 


where the primes indicate that the wings are pitching about 
and the moments are measured about the leading edge. 
These expressions may be inserted in equations (4) and, since 
the integrand becomes zero everywhere except at the point 
I, =t, the contributions of the initial pulse to the expressions 
for the lift and pitching-moment coefficient developed by an 
arbitrary variation of a and @ with time are 


Te? s 
may, tape? 
C. Pune 9 T0 3 
me AV. 04 V 


where œ and @ are evaluated att’. 
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FIGURE 2).—Undegenerate form of Initial pulse 
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The variation for t^ 7»0.— The integral equation (38) is still 
perfectly valid when applied to the flow field for /^70. 
It is convenient to rewrite the equation in this case, however, 
so that the effects of the vorticity on the wing and in the 
wake are separated. Thus, 


| Ant" (1) 
2-7 Je-Vgt^ T—34 


dz 
(44) 





e- bgt’ d c 
| ea f Ault, Ti) fy, 1 
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where Au* (qı) is the value of Av in the starting vortex wake. 
It is independent of /' since its value at all points along the 
line ab in figure 21 is the same as at the point a. 
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FIGURE 21.— Wing and trailing vortex sheet in zi^ plane, 


A reduction of equation (44) can be obtained for the case - 


of the sinking wing, where w= — Va, by using the inversion 
given by equation (39). Thus, 


rAu(", 2)y(z FPPC- VE A= A—4(22--2Y4'—c)Vwa-- 
f - ANu* (i) È Yin t Veret Ve) "T 


(z— 21) 


Since A is given by the relation 
s= (^ aut nd *  awriedd 
2 =[" Ault, 21) n= |", u"ridzi 


it follows that 


ya(2x +2 Fë —c) 1 
x (2+ Vat!)(e— Vet’ —x) 


og EE GET EHV dn, (48) 
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According to the Kutta condition Au(é’ x) remains finite as x 
approaches e— Vt; the integral equation for Au*fz,) thus 


becomes 
° ap® Lit V 


Vie= — 


which was derived and studied originally by Wagner (refer- 
ence 14). 

The section lift and pitching moment can be derived in 
terms of Au*(z) in the following manner. By definition, 


the section lift I is 
za FE Ag (4T 


(z+ Vet Ae Va — 2) 


41] 


Since the value of Ag is zero at the leading edge and at the 
trailing edge is equal to the total circulation T', two alterna- 
tive forms for the lift can be written 


ios c— Fot’ 
I= ps V + po i oe Ag ajde (48) 
and 2 
[= pe(c— vet) = » — Po 37 dt’. Vat’ Aul, x) dr (49) 


By substituting equation (45) into (49) and integrating, it — 


can be shown that 


* 
f= pocte [. — e 
F pot Fier- n c— Vot’ Ar Vo t’) (zx,—c4d- Vt) d 


bi | 


Since the value of Au*(ri) has been determined by Wagner, 


equation (50) can be evaluated to obtain a solution for the 
section lift. A plot of the section lift coefficient is shown in 
figure 22. Initially there is the pulse having an intensity 
defined by equations (43). 
value of the section lift coefficient starts at one-half its 
asymptotic value. It then increases, slowly approaching its 
asymptote of 2re. 
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FIGURE 22—Indicial Ift curve for incompressible fow. 


By definition, the pitching moment can be written 





— ies RAAPI TP da 
c— Vg : DA i 
-—n[. | (Fo! -2) SE dz (51) 


where the moment is taken about the leading edge and is 
positive when the trailing edge i is forced down. A develop- 
ment, similar to the one given for the lift, gives 


m-—2l | (52) 


After the pulse at //—0, the 
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This result, that for 4^0 the indicial center of pressure re- 
mains constant at the quarter chord throughout the motion, 
is classical. 

If the boundary condition for a pitching wing,.w,=— 
(z+ Vat’) 6, is substituted into equation (44) and the inversion 
given by equation (39) is again used, it can be shown in the 
same manner used in the derivation of equation (46) that for 
z-c—V,t^ the relation 


applies. This integral equation applies to a wing pitching 
about its leading edge. If, instead, the wing is pitching 
about the three-quarter-chord position, an essential simpli- 
fication is achieved. In this latter case, downwash is given 
by the expression 


Wy,=—é (2+ Vat c) (53) 


and the resulting integral equation becomes 


c zı-+ Vot 
o= f y,o at (zi) VEA, dn (54) 


where At4* (21) represents the vorticity in the wake following 
such a motion. The solution to equation (54) is simply 


Aus" (zi) =0 (55, 


From equation (55) it follows that the total indicial lift for 
t' ^0 on a wing pitching about the three-quarter chord point 
is zero, and that the wing wake is free of vorticity. Further, 
it can be shown that the total indicial pitching moment (still 
measured about the leading edge) is 





3 
m=— TE Vab | (56) 

The transfer of equations (55) and (56) back to the case 
in which the wing is pitching about its leading edge can be 
readily accomplished by means of the boundary condition 
shown in figure 23. Hence, if (e; DIA refers to the lift coeffi- 
cient on a wing pitehing about the three-quarter-chord point 
and (Cm ^)a4 refers to the pitching-moment coefficient meas- 
ured about the leading edge of a wing pitching about the 
three-quarter-chord point, then 


£ 
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(57) 
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FIGURE 28.— Change in boundary conditions corresponding to change !n pitching axis. 


By means of equations (52) and (57) the expressions of 
the three indicial functions, c,,’; Ci and Cm,» C&D all be 
written in terms of the indicial lift function fort >0. Hence, 


f 


Cn. = E e. 
à 
C, =] Ci, (58) 
c "TE: C EN 
me 16 la 8 


The variations of the four indicial functions will be shown 
later in Part IV. For values of ro(= Vot’/c) larger than those 
shown in figure 22 the approximate equation suggested in 
reference 16 can be used, namely, 


E „= 1- gi (59) 


2+ 


This alternative result has, according to reference 16, an 
error of 2 percent or less for the entire range of time froin 


0+ to infinity. 


INDICIAL FUNCTIONS FOR A TWO-DIMENSIONAL FLAT PLATE, Me=m0.5, 0.5 


When the Mach number is no longer small, the analysis 
in the preceding section must be modified. As an example 
of this modification, we shall evaluate the indicial response 
on a sinking wing flying at a Mach number equal to 0.5 and 
the indicial response on a sinking or pitching wing flying at a 
Mach number equal to 0.8. 

Since the wing is two-dimensional, the partial diff erential 
equation governing the flow field (equation (9)) reduces to 


Prt Par Pit (60) 


where it must be remembered that the axes are fixed with 
reference to the still air at infinity and the wing is moving 
in the 2—0 plane. The equation for the loading coefficient 
remains as in equation (12). The analogy which existed in 
the incompressible case between the theory for the unsteady, 
two-dimensional wing and slender-wing theory exists in this 
case between the theory for the unsteady two-dimensional 
wing and the theory for a steady-state, three-dimensional 
wing traveling at a supersonic speed. Thus, in the three- 
dimensional, steady-state case the partial differential equa- 
tion governing the flow is 


yy F 911 = 8^ e (61) 


and the equation for the loading coefficient is 


The boundary conditions are in both cases that o, is given 
over a portion of the z=0 plane. It is evident by a compari- 
son of equations (60) and (61) and equations (12) and (62) — 
that results from the three-dimensional, supersonic, steady- 
state case (hereinafter referred to as the steady-state case) 
can be transferred to the two-dimensional, unsteady case 
(hereinafter referred to as the unsteady case) simply by 
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replacing r, y, and B in the former case by f, z, and 1, re- 
spectively, and by dividing the result for the loading coefficient 
bv Mo. 

The analog to the boundary condition for the problem of 
finding the indicial loading on a two-dimensional wing fying 
at a subsonic Mach number (fig. 24(a)) is the boundary 
condition for the problem of finding the loading on a constant- 
ehord, swept-forward wing tip with a subsonic trailing edge 
such as that shown in figure 24(b). The Mach cones in the 
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iat Unsteady case. 
(b) Steady-state case. 


FiGr &E 24, — Boundary conditfons for a two-dimensional unsteady wing mov- 
ing at subsonic speed and the analogous three-dimensional, steady-state 
wing. 
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steady-state case, traces of which are shown as dotted lines 


in figure 24(b), become, in the unsteady-state analog, the E 


locus of the sound waves which started at t=0 from the 


leading and trailing edges of the two-dimensional wing _ 


(fig. 24(8)). Finally, the analog in the steady-state field _ 
of the unsteady wing would be a flat plate for the unsteady 
sinking wing and a plate with a linear variation of twist for. 
the unsteady pitching wing. 

Just as in the section on incompressible flow, the analysis 
wil be divided into two parts. In eases for which 14,750, 
however, the indicial functions contain no pulse at £—0. 
Hence, the first part of the study will be concerned with 
the behavior of the indicial functions in an interval for 


which ¢ is small but finite and the second part, with their B 


asymptotic behavior. | 

The early stage.— The analog which exists between the 
steady-state and unsteady cases may be utilized to great 
advantage since the special methods and techniques developed . 
for the solutiou of problems in the former case may be applied 
to the solution of the analogous problems in the latter 
field. In this manner an exact solution for the loading over 


the first five regions shown in figure 25 was obtained fora — 


Mach number equal to 0.8, and for all the regions indicated 
for & Mach number equal to 0.5. Solutions for larger 
values of the time could also be calculated, but the labor | 
involved in computing such cases becomes prohibitive and, 
as will be shown later, approximate methods can be developed 
which extend the solutions for the indicial lift and pitching- 
moment curves to their asymptotic values. 
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FiGrRE 25.—Reglons used [n analysis of subsonic unsteaty wing. 


The analysis used to calculate the loading over the regions 
shown in figure 25 is outlined in appendix A. Plots of the: 
indicial loading on a sinking and a pitching plate flying at a 
Mach number equal to 0.8 are shown in figure 26. At 
time equal to zero the loading is constant for the sinking 
wing, and as time increases the loading-coefficient curve 
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(b) Wing pitching about Its leading edge. 


FIGURE 26.—Varlation of two-dimensional indicial load distribution with percent chord 
for à Maoh number equal to 0.8. 
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approaches the familiar two-dimensional, steady-state shape 
given by the equation 


(63) 


Ap, de [ome 
Qo B I 





where z is distance from the leading edge which is at r=u._ | 
Figure 27 shows the load distributions for Ma=0.5 and 0.8 
at the last value of time for which the exact loading curve 
was calculated. Notice that in each case the distribution is 
essentially the same as that obtained at time equal to infinity 
(i. e., the agreement is good with the curve produced hy 
multiplying the right side of equation (63) by a constant 
` factor}! "The use of this fact simplifies the subsequent 
analysis concerning the asymptotic behavior of the indicial 
curve. FP sj l 
The indicial lift and pitching-moment functions were also, 
calculated (see appendix A) up to the time 75—2.333 for the — 
wing flying at Ma=0.5 and ro=4 for the wing flying at 
. My=0.8. Their variation in this interval will be shown in a 
subsequent figure (fig. 59). It is evident from a glance at 
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(a) AfemO.5. ; 
(b) Afom0.8. . 


FiGvRE 27.—Load distribution at end of early stage. 


V A similar result was noted in the study o {the load distribution on swept-back wings with subsoaic leading edges (reference 17). 
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the figure that the calculations must be extended beyond 
this early stage since the asymptotic values are not even 
closely approached. 

Before studying the nature of these curves for large values 
of time, however, it is useful to examine them with reference 
to the discussion in the previous section on incompressible 
flow. For example, it was’ pointed out that the indicial 
center of pressure on the sinking wing remained at the 
«quarter-chord point for time greater than zero. Hence, 
let us eonsider the location of the center of pressure on the 
sinking wing when the Mach number lies betw een 0 and 1. 
By means of the indicial curves for c; and Cm, ' and by the 
relationship 


Cm, = iS (z/ C)e, pt, 


the variation of (x/c),,. is easily evaluated ((z/c),,. is the 
distance between the leading edge and the center of pressure 
divided by the total wing chord). This variation is shown 
for the two Mach numbers in figure 28. It is apparent that 
the eenter of pressure is very close to the quarter-chord 
position for values of time greater than those for which the 
exact calculations were carried out. In other words the 
significant effect of compressibility on the location of the 
center of pressure is limited to the interval 0<%,<2 for 
jf,=0.5 and to the interval 0<é,<5 for 1/,=0.8. 

Likewise, it is apparent from the discussion of the incom- 
pressible case that the indicial functions for the pitching 
wing can also be expressed in a more convenient form by 
shifting the axis of rotation from the leading edge to the 
three-quarter-chord point. The values of (c;),4 and 
(Cm, ‘Jaa for 14,—0.8 were calculated from the definitions 
given in equations (57) and are shown in figure 29. In- 
spection of this figure again shows that at a Mach number 
equal to 0.8 the compressibility effects are limited to the 
interval 0«£&« 5. . 





g 2 4 6 8 
fo 
Figt 22 28.—Center-of-pressure variation on sinking wing during early stage. 


The later stage.—It follows from the preceding discussion 
that when č is large, the values of the indicial functions 
Cn, ^. ¢i,’, and Cm’ for compressible flow can be expressed in 
terms of c, by equations similar to equations (58) which 


were derived for incompressible flow. Thus, after several 
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FiGU RE 29.— Variation of lift and moment on pitching wing during early stage; Mew 0.8 


chord lengths have been traveled, one can write, on the 
basis of the asymptotic values shown in figures 28 and 29, 


Cx, = — C1 /4 
C, —36i [4 
"== —(3e,,/16)—(x/88) 


(64) 


Cu, 


[d 


It remains, therefore, to determine the asymptotic behavior ` 


of e. - 

Consider the steady-state solution for the lift on & two- 
dimensional, flat lifting surface traveling at a subsonic 
Mach number. As was pointed out by Wieghardt (reference 
18), if the lift on such a surface is represented by placing at 
the quarter-chord point a vortex which has the same circu- 
lation as that developed by the wing, the angle of attack 
measured at the three-quarter-chord point will be the same 
as that of the flat plate. Extending this concept.to include 
the unsteady effects, an investigation will be made of the 
variation with time of the vortex strength which will main- 
tain a constant angle of attack at the three-quarter-chord 
station following an impulsive start at ¢’=0. 


The analogous problem in steady-state theory becomes one m | 
of finding the strength of the vortex system, shown in figure © 


30, which gives a constant value of w along the line CD. 
In the vicinity of the origin, of course, this representation 
gives & poor approximation to the original boundary-value 


problem. On the other hand, in this vicinity the exact values — 


of lift and moment have already been determined. 
Esch vortex composing this system lies along the line AB, 
extending from minus infinity toward the origin, and trails 
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back parallel to the xı axis to form the trailing vortex sheet. 
Note that, for convenience, the origin of the axis system hes 
been located at the quarter-chord point. The solution to 
such a problem in steady-state, lifting-line theory would re- 
sult from the solution of the integral equation ™ (for a de- 
velopment using the notation adopted here, see reference 19): 


Be Agy) dz, 2 
re | an [(z — z} — gXy — y) 


fe 


vx, ez - Mot, 
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ye gee 


í For fe Ill heef 
VAR 


(à) Unsteady ease. 
(b) Bteady-state analog. 


FtGuRE 30.—Vortex systems for two-dimensional unsteady and analogous 
three-dimensfonal, steady-state wings. 


i4 T'he symbols f and T gro used to Indicate that the finite part Is to be taken. 
= f) 

z. 
. Yi-y 
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Thus (see reference 18 or 20), T a-p -E 
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where Ag is not a function of z; since the strength of a trail- 
ing vortex is, of course, constant. "The area of integration 
r is the region within the forecone springing from the point 
x,y. 

If the above equation is transformed by means of the 
analogy to represent the solution of the unsteady problem 
(sco fig. 30(8)), B, z, and y are replaced by 1, é, and z, rc- 
spectively ; and Ay, the total jump in pátentinl at a given 
section, is replaced by the circulation T. Hence, 


7-30 gos ye zia 7 


where r, as indicated in figure 30 (a), is the area in the forc- 
cone from the point P which lies always along the line 
z—(1/2)0— Aft. Integration with respect to & reduces the 
last equation to 


e A —— 
ETIE AE FrN o. _., 
M Ts —! 253) —(z—zjy 





which, by means of the substitution 2,/e=—2,/2 becomes 


along the line r= 5— My 


lu. 0 Otea Me) 00 
UI EB Jo (ua— Ào- 23)4 ((Xo— za) Aot ui— T3) 


where Ao = 2A Loto— uo, uo Mol (1 +A), = 2AL,/(1 —My), 
p=l/(1-HMo), &y—t/e, and where, of course, w is & constant 
equal to — Foa. 
A solution for (aa) in the integral equation (65) may be _ 
obtained by expanding T in a series of the form 
l | (66) 


is aem Vds CES MUR os EE 


Place equation (66) into (65) and expand in powers of 1/X,. 
There results the expression 


drz (6 5). 


w cmn L7] . (R7 
E -rat MM Cx. UMS A ee 
in which | 
ao, 1 
ahy 


Hence if ao and 5, are chosen so that e; is zero, an expression 
for T will be obtained which represents the solution to the 
integral equation (65) correct to the first order in I/f (i. c., 
1/Aj) for large values of &. Further, if equation (00) is 
expanded in powers of 2a, there results 


mye j- es [4 (^-2 +... | 


L* f(y) dy ra en = G(r,ó)—G(z,a) where G (x,y) Is the Indefi- 
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which becomes, using the condition for c; and relating z, and 
h by the equation of the leading edge, 72=2Af4—1= 
(2Uu'/c) —1, i 
| aCaV, = 1 
rg |: zm | (68) 


By choosing tlie values of a, a, and 6; given as follows: 





5; =(ao/2)— (1/8?) 
0.5 50 28 26.67 (69) 
S 36 27 15.22 


and placing the resulting expression for T'(z;) as given by 
equation (66) into equation (65), the values of —1w/Voe 
shown in figure 31 were obtained. This figure demonstrates 
the accuracy to which the first two terms. of the series 
expansion for I'(a3) yields a constant value of w for the 
constants given in equation (69). — 

The relation between circulation and lift has been derived 
and presented as equation (48). This expression can be 
written : 


1 
l= poW oI + pot i [Sero Xod Lo 


where x—z/c. In order to obtain a complete expression for 
the section lift, it is necessary to know the chordwise varia- 
tion of Ag. Since equation (66) gives only the total vorticity 
and not its chordwise distribution, some assumption as to 
' such distribution must be made. In lieu of this, the result 
presented in figure 27 suggests that for large values of time 
the value of Avtro,a%) used in the equation for section lift 
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(a) Momy.5. 
(b) Mamü.8. 


FIGURE 31.—Variation of downwash at three-quarter-chord position. 
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FIGURE 32.—Chordwise variation of circulation at end of eariy stage; 3/a=0.8. 


can be expressed by the product of Ag( œ te) and T'(r;)/T ( e). 
In other words, for large values of ra the shape of the chord- | 
wise distribution of vorticity is the same as the two-dimen- 
sional, steady-state value. An indication of the accuracy of 
such an approximation is shown in figure 32 for a Mach 
number equal to 0.8 where the precise value of Ag(rp,20) is 
compared with the approximation at r9,—4. 

Since T'( o) zal eB, the substitution of 


Agro, Fo) — Ae( o , zo)T(ro)/T( e) 
gives for the section lift coefficient 


TI ,3 1 dT = 
6—yc 9 Pe dt (70) 
By means of equations (66) and (70) and the values of the... 
constants given in (69), an expression for c;, can be written 
which is valid for large values of ro This expression is 
somewhat cumbersome, however, and it is difficult to apply 
in subsequent analysis. Hence, it has been replaced by a 
simpler equation that is equivalent " to three decimal places 
for values of ry greater than about 10 for both Mach num- 
bers 0.5 and 0.8. These equations are given by the following: 


for Af,p=0.5 | 
. 4 _ 44.218 i 71) 
: 35-F2n (5-25 


2T 4 l 
———11 
« p 


and for .4,20.8 !'5 


Qe 1.736 


70.83 - 
"« B usa a | (72) 


Ci 


IT In order to facilitate fairing the Initia] portion of the indicial lift curve mto the approximste solution obtained for the subsequent variation, the latter results were shifted slightly in terms 


of Vat'[c. This did not affect the asymptotic behavior of the curve. 


i5 The value for er given by equation (72) disagrees slightly with the value given In the superseded TN 2103 (see footnote 1}. The methods for caleniating the results are somewhat —— 


different, and the value given here is considered to be more accurate. 


2124832—354— —28 
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The equations for e,,’, e and M for r; 10 can be 


calculated from these expressions by means of the relations 
given in equations (64). Equations (71) and (72) will be 
considered valid for values of v; greater than 10 for both 
Mach numbers. Values of the indicial functions between 
t=10 and the highest 7) calculated by the exact method 
were approximated by interpolating between the two results. 
Tabular values of all the indicia] functions caleulated for 
the interval O< 74,10 are given in tebles I and II and 
curves are given in Part IV. 


TABLE I.—INDICIAL LIFT AND MOMENT FOR TWO- 
DIMENSIONAL SINKING WING. AM$0.5 
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TABLE II.—INDICIAL LIFT AND MOMENT FOR TWO- 
DIMENSIONAL SINKING OR PITCHING WING. 21,—0.8 





Vot'jc Be, /?* Aem y [x 
0.0 0. 478 0. 955 0.318 : 
.l . 466 om. - 814 - 61 
.2 . 454 .880 31S .. 50b 
3 .442. . 830 . A106. - .D81 
T . 430 -730 . 328 . 57b 
-6 . 423 . 680 .à39 . 582 
.6 . 426 - 05d . 955 . 595 
7 433 . 632 M P" . 608 
.8 .. 442 .020 ; . Bb .821 
.9 451 .815 .897 . 038 
1.0 . 461 - G11 411 . 644 
1.5 . 907 .612 PE YE: . 688 
2.0 . 546 .019 . 525 -TIS 
2.5. . 581 . 628 . 566 . 149 
3.0 610 .041 . 608 » (3f 
8.5 . 032 . 655 „~ 625 . 746 
4.0 . 652 . 670 . 648 ^ .755 
4.5 . 670 . 684 „005 . 7 
5.0 687 . 693 . .881 . 778 
8.0 714 . 121 . 710 . 190 
1.0 .798 — . 743 . 736 . 806 
8.0 -T . 763 . 160 .B21 
9.0 . 770 . T82 .T8D -. 835 
10.0 . 708 . (98 198 . 849 
e 1.000 1. 000 Loo 1. 000 
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INDICIAL FUNCTIONS FOR A TWO-DIMENSIONAL FLAT PLATE, Mi=1,0 


The general results, obtained in the preceding section for 
the early stages of the motion and presented in appendix A, 
for the indicial loading over the sinking and pitching wing 
may be extended to the sonic case. Furthermore, the two 
intervals for which analytic results in a closed form were 
presented in appendix A now cover the complete time range 
since o< ro S Mif (Y -- My) becomes 0< ro SS 0.5 and My(14 
My € oS Mo/(1—M,) becomes 0.5 rj o. Hence, by an 
appropriate limiting process, equations (A8), (A9), (A10), 
and (A11) become for 0€ 74, 0.5 . 








Ci, = 
OLENS Oa =EN 
l l | (73a) 
€, 2d ro 
me == (4/8)— (2/8) re 
and for 0.65 mS 9 — 
ti s) (22 ro—1+are cos —2— =) 
Cm '=-@/ 7) aki V2To— i+ are cos | 
0 
M p [re To (124 arc cos at] (73b) 
0 = 
„=ar ant I2tg9—1+ 
1 (1 n) arc cos t~. -] 
a a ^ To | 


Since the magnitude of the functions in equations (73b) 
grows indefinitely with increasing time, the assumptions of 
linear theory are eventually violated. However, for moder- 


“ate values of rẹ these functions have the same order of 
magnitude as similar indicial curves for Mach numbers other 


than 1. "These effects are illustrated in Part IV. 


n FUNCTIONS FOR A TO med QE PART PLATE, Him 1a 


-- "7 4 -A- 


The method of obtaining solutions for the atic functions 
at supersonic Mach numbers parallels the development 
presented for the subsonic Mach numbers. The steady- 
state analog to the supersonic unsteady wing problem is a 
consíant-chord wing tip with a supersonic trailing edge. 
(See figs. 33 (a) and 33 (b).) It is well known that the 
problem of finding the loading over wing plan forms with 
all supersonic edges is one of the simplest in three-dimensional, 
lifting-surface theory. In fact, since the upper and lower 


-— T— LL 


surfaces.&re noninteracting, the solution is determined by _ 


integrating sources within the Mach forecone. The analysis 


for c;, has already been carried out in reference 21. 


The analysis used to calculate the loading in terms of 

zo —z/e and ty=t/e over the three regions shown in figure 34 is 
outlined in appendix B. Anexample of the manner in which 
the loading varies with time over a sinking or pitching wing 
traveling at a Mach number equal to 1.2 is given in figure 35. 
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(b) 


(4) Unsteady case. 
(b) Steady-state case. 


FrGURE 33.—Boundary conditfons for a two-dimensional unsteady wing moving 
at Supersonic speed and the analogous three-dimensional, steady-state wing. 


The expressions for the mdicial lift and pitching-moment 
coefficients are given analytically in appendix B, and plotted 
in Part IV. Tt can be shown that the results given in appen- 
dix B reduce to the expressions given by equations (73) 
when Jo is allowed to approach 1, so that there is no dis- 
continuity in the theory in passing through the sonic range. 


SUPERSONIC STEADY-STATE LIFT 


Since the next section contains & problem involving a 
complicated acoustie plan form, it may be helpful to consider 
first a problem involving a very simple acoustic plan form 
but otherwise similar to the subsequent analysis. Hence let 
us inspect, using Kirchhoff's formule, the problem of finding 
the steady-state loading on a two-dimensional flat plate 
traveling at & constant speed (equation (15)). 


Since the upper and lower surfaces are noninteracting, - 


we ean use the special form of Kirchhoff's formula given as 
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Fict'2& 34.— Regions used in analysis of supersonic unsteady wing. 


equation (28). In the plane of the wing this equation be- 


comes 
1 LL W.drddy, 0 
e] Se J v(z—2)?+(y—y)" 
Figure 36 shows the positions of the wing in the zy plane. 
The wing has constant speed for £50. The point P(z,y) 
is chosen on the wing and ahead of the wave which started 
at time zero; therefore, P(z,y) lies in the region which has 
attained its steady-state value. Further, the value of w, 
is constant over the acoustic plan form in such a region. 
This constancy reduces the problem to one of integrating 
[r—nY-F(g—45Y] ^ over the ellipse representing the 
acoustic plan form. 
The equation for S, can be determined from equations 
(20) and (21). In this case, equation (21) becomes simply 


(74) 


tı = — Ahr 


Eliminate r between this equation and equation (20) and 
there results in the z=0 plane 
2 2 z, V = 
(=r +y yd - (H1 75) 
AL 
That equation (75) Is the equation of an ellipse with one focal 
point at z,y can be readily verified. It is more convenient, 
however, to change to a polar coordinate system with origin 
at P. Henceset l 


r—t=r cos @ 

y—y =r sin f 

dy, dz, —rdrdó | 
Then equation (75) becomes 


Mor = Mtt r—r cos 6 
or 
` r4 Mg 


ipn) +cos 8 (76) 
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Ti ae 
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and, inim, equation (74) becomes simply 
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The integral is not difficult to evaluate So 
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and, finally, by equation (12) 


oa Ap _ 4 = uM, )-3m* = uS 
q FoMo yMê— 1 vM—1 | 


which is the familiar Ackeret value for the loading on a- 
two-dimensional flat plate. The lift coefficic:.t, of course, 
follows immediately as 


. a| ^ 





4a | 
o ROT " = -- (78) 


INDICIAL LOADING FOR SINKING TRIANGULAR WING WITH SUPERSONIC 
. EDGES 


The seven dedins Tie analytic expression for the indi- 


each of seven regions. These regions are determined by 
the positions of th» various wave fronts relative to the wing 
plan form (fig. 37). For t«C0 the wing is motionless, its 
leading edge lying along lines represented by the dashed 
lines in. figure 37. The trailing edge is considered to be 
supersonic and hence its position is immaterial, since the 
solution can be cut off wherever desired. At f=0 the wing 
starts suddenly to move, and for £750, travels forward at a 
constant speed Vo. After a certain time ¢ has elapsed, the 
wing has traveled to a new. position, also shown in the figure. 





Percent chord 


(a) Sinking wing. 
(b) Wing pitching about its leading edge. 


FicvRE 35.—Varlation of two-dimensional indicia! load distribution with percent 
chord for a Mach number equal to 1.2. 





cial loading over the triangular wing has a different form in *.—. o 
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In this same interval of time, pressure impulses have trav- 
eled out in spherical waves from every point of the region 
whieh the wing has occupied. "The trace on the wing of the 
sphere starting from the wing apex at (—0 forms the exter- 
nal boundary of region 7. "The area outside this circle and 
within the traces of the cylindrical waves (the envelopes of 
the spherical waves) generated by the leading edges at £—0 
forms region 4. Region 5 is formed by the overlappping of 
these cylindrical waves, and the solution for loading within 
it can be found by a suitable superposition of the solutions 
for regions 3 and 4. Region 1 lies between the cylinder 
trace on the wing and the leading-edge position at time i; 
the loading in this region cannot be affected by the manner 
in which the wing started its motion since it lies outside the 
starting cylindrical waves. Hence, the loading in region 1 
is the same as that on a swept wing flying at a steady super- 
sonie speed. The solution in region 2 can also be obtained 
from steady-state lifting-surface theory, but, whereas in re- 
gion 1 the feld is two-dimensional (i. e., invariant with dis- 
tance measured parallel to the leading edge), in region 2 
the field is conical. Region 6 is formed by the overlapping 
of regions 2 and 4. Finally, region 3 is that area com- 
pletely unaffected by waves from the wing edges. In the 
following subdivisions the analysis of each of the separate 
regions will be discussed. 

Region 1: The loading in region 1 of figure 37 is equal to 
the loading on a two-dimensional flat plate moving at a 
constant velocity given by the component of stream velocity 
normal to the leading edge of the triangular wing. Since 
this component is supersonic, the loading is of the Ackeret 
type and is given by 
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FIGCRE 3¢.—The seven regions used in the analysis of the triangular wing with supersonie 
leading edges. 
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QD “Fat (79) 


where 8 —4/ Af,?— 

Region 2: The A loading on & triangular wing 
with supersonic edges has been given by several authors (see, 
for convenience, reference 22) so the expression for the load- 
ing in region 2 can be written immediately for the coordinate 
system shown in figure 38 as 


B'my—(rrAMg E 
aaa tere sin Bima + ALS — vi 


s TV a i 
EA 
are Sn aima FMA Ay EO 
Region 3: Since region 3 is unaffected by the edges of the 
wing the solution for the loading therein can be written as 
in reference 21 
£3! E 


Region 4: The solution for loading in region 4 can be 
obtained from consideration of & two-dimensional wing 
starting from rest and moving with velocity V, normal to 
its leading edge. "This problem has been treated in refer- 
ence 21 and the solution written there cau be written for the 


right-hand side of figure 37 as 
VAL 2-1 
OM. (F+are amp 


Ao, | 
= arc cos 
where the notation, as defined by figure 39, is 


(81) 
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FIGURE 33.—The rave coordinate system. 








and since 
cin A=m, sin A= l 
41 ix 
cos A= E NER x: 
414m? 
Then 
ma-— 
PEEL g : 
414m 


The equation for loading now becomes, in the coordinate 
system of figure 38, 


(2) | are af m M ymz-—|yl--t1 +m) 
gems Jy1-4-m mz |y R Me) 


yB*m*—1 icu arc sin au (82) 


"^M, 0 

Region 5: The solution for loading in region 5 can be 
obtained by superposition of the solutions for regions 3 and 
4. (See fig. (40).) If the solutions for the two sides of 
region 4 (obtained from equation (82)) are added, the result 
gives twice the required value of w, on the wing, as well as 
undesirable pressures off the wing. However, subtraction 
from this sum of the solution for region 3 (equation (81)) 
reduces the downwash w, to the proper value, and also 
cancels the excess pressures. ‘The resulting expression can 
be written : 





| mM (me—y)+O1 tm} 








x 5 eS » Em (mr —y +m My) 
m M (me+ty)+a+myt 
ido Vi --m*(mz Ly 4 m M 4f) 
EA 1 MITY mz—y 
P ( aresin TEE, "Nr S aresin T J| (83) 


Region 6: The loading in region 6 can also be calculated 


by superposition. To find the loading in this case, add the 


REPORT 1077—NATIONAL ADVISORY COMMITTEE FOR, AERONAUTICS 


— Figure 40.—Solutions superimposed to obtaIn equation (83). 


solutions for regions 2 and 4 (equations (80) and (82)) and 
subtract the solution for region 1 (equation (79)). There 
results 
G Ap) _ — fat — fares n É'my— (+M 
T in Smi F M,)— y 
B'my +(x + M) 4 
B[m(z 4- M o£) +y] 


mE mMy(mz — |y[) 3-1 -- m? 
À TUE (mx—|y|-+-m M Vt) Em. 

| | ymi .. ma—lyl\} 
mM, (Frere sin 5 


arc sin 


(84) 
Region 7: The solution for the loading in region 7 can be 


obtained by means of equation. (28). The analysis used in 
finding the solution in this region is not diffieult but the. 


- algebra is rather involved. It is useful at this point to intro- ] 


duce polar coordinates (see fig. 41) such that 
z—1,—r cos 0 
y—y =r sin à (85) 
da dy, —rdr d 


From equation (85), equation (28) can be written in the 


form 
| ope, farao 


The acoustic plan form for points in region 7 is the region 
bounded by three curves as indicated in figure 42. The are 
between @, and 6, is determined by eliminating T between the 


equations 
| r-(— T) 


. (80) 
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(the equation for the inverse sound waves) and the equation 
for the left leading edge 


y= mat} 


The arc between 0, and 6, is found by determining the acoustic 
intersection (eliminating T) of the right leading edge with 
the inverse sound ware; and, finally, the arc between @ and 6; 
Is given by the equation r—£. The equations of these arcs 
can be written in polar coordinates as 


LU. 
mcos 0--m M,-Esin 6 Ses 8 


rz 04 «0 «0, 


|. m(r--Md)y . lo E8x, 
. m cos 6-+mAf,—sin @’ 








S (2, Y) 
Jy, 
Zi 
FIGURE 41.—The r$ coordinate system. 
——— Wing at to 
Wing now 
8; 70, -2x 
Yi 
/ 
Á iN 
Á M i 


X 


FIGURE 43.—Acoustic plan form for point In region 7. 


Using these expressions, equation (86) reduces to =i 


Voa ^s m(t +M) Fy d6-- = 222 
2r Je, m cos 6+mAs,+ sin 0 


Vou " td6-I- a 





2v Je 





Via f m(x--My)—y 25 
2r Je, m cos 09d3-m M,— sin 0 


and taking the partial derivative with respect to £ (to deter 
mine the loading according to equation (12)) one finds ?? 


Ap\ 2em (^ — do (a (^ 
cm a jo, mAf,+m cos 6+ sin trl, d64- 
2am (A dé 
rz m, +m cos 8— sin @ (87) 
In evaluating this equation, the following integral is used: 
for —7S8€X- E 


dé | 2 
J mi,]-mcos8-Lsin8 /gimi—] 


m(M,— 1) tan (8/2) -- 1 


Jami "" 


arc tan 


Since equation (88) is valid only in the interval—x €6& r, 
care must be exercised in applying it because the angle 6, | 
may be greater than r (as in fig. 42). In case r S6, it is 
convenient to introduce the angle 0,'—0,—2z. "The expres- 
sion for Ap/g, can then be written in two forms, according 
to whether 6, is less than or greater than r: 


for yz0, 4, « 
P) 4am m(Jfa— 1) tan (6,/2)—1 i 
SE yi tin Le ee MEL e 
t 7T mB m’—] | = tare us 4B*m*— 1 


m(Mq— 1) tan (6/2) 1 | 
Jami 
m(Af,— 1) tan (0/2) +1 — 
Jem? I 
m(Af,—1) tan ga: 4 2% 6) (Boa) 


arc tan — ———É—————— rif, 


Jemi 


are tan 


arc tan 


for y> 0, 7 € 6, 





dp) tem _ [ore tan Mat) tan (6/2)—1 _ 
o/z xy f?m’—l] | are i J8*m1— 1 
m(lfo— 1) tan (6/2) 1 pes 


arc tan — + 
JB'mi—1 


m(M,—1)tan(0/2)--1 _ 

o Pma 

m(My— 1) tan (8/2) —1 
VBimi—1 


arc tan 


2 
bk war, 5 2x) 
(89b) 


arc tan 


*' The limits f, 61, and 6; are all fanct!ors of ¢ but in moving the partial derivative through the Integral sign the terms Involufng h/t, 0630, and 06 X all cancel one another, 
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where 
— unu Lego aed 
(0« 0, v) 
6,—arc cos nutno- t m Eu 
(0«C6,« v) 
masse hU ne AERE 


9, = 0,— 2f T < 04) 


The limitation on 6y can be given both an analytic and 
geometric interpretation. Thus, equation (898) applies for 
O<m(rti)<Sy and equation (89b) applies for 0€ ym 
(z--t). These regions are shown in figure 43. Because of 
the geometrical symmetry about the xz axis, equations (89) 
suffice for the determination of loading throughout region 7. 

An isometric drawing of the load distribution on the right 
panel of a triangular wing with supersonic edges is shown in 
figure 44. The positions of the spanwise sections were 
chosen so that each of the regions 1 through 7 is represented. 
It is to be noted that the results for region 7 show no unusual 
characteristics and, in general, the distribution is similar to 
the steady-state loading on a triangular wing. 





r- Equation (89(8)) 
! valid 





DOCE 
unos 





~ -Equo tion 
(G9(b)) valid 


—--—-Wing af tsa 


Wing now 





Fiol RE 43.—Regions in which equations (89a) and (89b) are vaiíd. 


INDICIAL LIFT AND PITCHING MOMENT ON A TRIANGULAR WING WITH 
SUPERSONIC EDGES 

The indicial lift and moment could be obtained, of course, 
by integrating load distributions caleulated by the method 
presented above. However, it is far simpler to use the 
methods outlined m the section entitled *Boundary-Value 


Ap 
Gok 


Fiaurg 44.— Distribution of {ndicial loading on right pane) of supersonic-edged, iri- 
angular wing. _ 


* 











—— — — froces of regions shown 
in figure 3 


Problems Involving Noninteracting Surfaces" in Part II. 
In particular the second of the two methods outlined therein 
will be applied. 

It was seen in the above-mentioned section that the lift 
and moment coefficients for a supersonic-edged triangular 
wing could be found by solving a related steady-state lifting- 
surface problem. For the cases of indicial sinking and pitch- 
ing, the boundary conditions for this related problem can 
be found readily from equation (34). Since y, and y, 
become, for z=0, the right and left leading edges, respectively 
the boundary conditions are given by 


E z = — 2? Vam (z + My) (90) 
&/s-0 
for the sinking wing, and 
ac) anode (z + Mg? (91) 
PA z= 


for the pitching wing. In addition, it was shown that the 
lift and moment coefficients for the wing in unsteady motion 
are obtainable by integration of the quantity 0@/dt. For 
convenience, this quantity is expressed as a chordwise load- 
ing factor in terms of the following notation: 
m AL, 3 Ap à í " 
Bem x dy (sinking wing) (92a) 
MoV, (tA v 
peti =P dy (pitching wing) (92b) 
. 2876 Js go 
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where 


s= mix +A) (93) 
It is found that (equation (35)) 


im Ó * Gri M) dxdt, 


y= ——— 
^os Jed yiia , 


(94) 


im? Ó 


(zxi- Md)? d xdi 
Ursi ot. 


LL ————————— (95) 
y —&y — (2) 
The lift and moment of sinking and pitehing wings will now 
be obtained by the use of P, and P\. 

Wing with constant angle of attack (sinking wing).—The 
solution for the chord loading on a flat, supersonic-edged, 
triangular wing starting from rest at (—0 and flying at a 
constant speed and angle of attack is given by equation (94). 
With the transformations 











J—J[p—24. 
i—i=& 
this becomes 
: x 
no Mg; 
P= m B [4t dr. 
T3 at Te tre 


This integral can be evaluated and gives for x<é (region A) 








P4 (96a) 
for —i<2St (region B) 
P= MS ree cos ( — 3 are cos e (06b) 
and for r<—é (region C) 
Pis (96c) 
where the regions are shown in figure 13. As has been 


pointed out, equations (96) can also be obtained by inte- 
grating the equations for the loading given in the preceding 
section. These integrations were carried out (in some regions 
numerieally) and the results were found to agree with those 
of the present analysis. 

It js now possible to write the indicial functions Cs, and 


C m in in the form 


x 


9 tg — M yl 
Cm gr, N ví DG AL Pad (97) 


c} Cg — Af gt 
( a JE ES OE BERE 
a SeqM, Ls 


where c, is the root chord, S is the wing area (equal to meg?) 
and the prime indicates that the pitching moment is measured 
about the apex, the positive moment being one which causes 
the trailing edge to sink relative to the apex. 


mi -- May Pur (08) 


Inter val 





FiGrRE i5. — Time Intervals use In expressing Hft and moment on unsteady, saper- 
d, triangular wing. 


Combining equations (96) and (97) one finds for the first 
interval shown in figure 45 


2 | -t 

CL. 3T oco! ME 
j E Pyr? 

s: NES Ez REDE cos (b 


Msg ve 
B arc COS 1 Jf dr H 


£M, 





(r+ M)dz4- Tm 





FEN l 
4(z d Md | 


This -— integrates to give, if fy=¢/eo, 


for OSASI I1 (first interval) 





—* (iple 
Cr = OH tt) (99a) 
E : 
Similarly for 4 TA = 1 ——T ty Ew S (second interval) O 
-- ls Mgv—1 | EN 
NES A (1 ts fo ) are cos 77 + | 
Ze are cos (.4.f,—t87)-+ 
zB 
o Mdo ips 70 — Arg (99b) 
rM, 
and for 1 T- 1 <t (third interval) 
C = e d E (990) 
R 


In the same manner the values for Ca io the various 


intervals can be determined by combining equations (96) 
and (98). There results 


for 0 €to X 32—— Af, + 1 (first interval) 


fn Bf pg is | 
Cn.’ = gaz (1 +5 Me) (1002) 


MT (second interval) 


8 


On, — Mr | 18 B—Mety—Mete!— 2h Vig MAF 


Me 6%) | 


(100b) 


5 2 HMot® arc cos sald ot TE D ~ arc cos (M.,— 


38 


and for M. T <t, (third interval) 


(100c) 


Numerical results will be presented in Part IV. 

Wing with linear angle-of-attack variation (pitching 
wing).—The solution for the average load on. a flat, super- 
sonic-edged, triangular wing flying at a constant speed and 
pitching at a uniform rate @ about its apex is given by 
equation (95). With the transformation 

L— t= 
i—i =h 
this equation becomes 


_ 4m? d f pee) Mss) 
P= P ot JGI— rr 
and the evaluation of this gives (for the intervals defined in 


fig. 13) 
for z&€ t (region A) 


n | 1+3 eaa) | 


for —£&z €t (region B) 


-i C (rM) 
2 (rM 


M, CHEM, 13:42Mj 
queen Sat GI Mpe VE | aor) 


di,dz, 


(1018) 


arc cos (-5) + 





and for x< —t (region C) 


P= B (101c) 
The equations for the indicial functions Cu and C, , can 


be obtained from the equations 


CU t 2 fom Mot : 

ra Ou See auy POTM Pde (102 
9) 

Us / 


= 2 amore a 
PUN RIT A ll m(z-J-MqgyjP.dz (103) 
S 
where the primes indicate the wing is pitching about and the 
moments are measured about the leading edge. 


_lytical form of the loading equation is different. 
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A combination of equations (101) and (102) gives for the 
lift nae 














for 0 StS i, - 1 (first areis 
C cupid '- Md) (1042) 
Ba 3 M, 2 : 
TN 
for MSS z 7 eee interval) 
CL =< I tes Mg + 3) arc cos Me ~- 
My, 
38 ? are cos (M,— 8t) + 
(2 22th ut ne Mg +41 4 Mg y | (104b) 
g° "gue Ig Pp re TNT 
` and for 7 o «t (third interval) 
C 28/38 (104e) 


Similarly a combination of equations (101) and (103) yields, 
for the pitehing moment about the apex, the results 


for 0 <ta € ——— Af ii (first interval) 


(1052) 








EX ——MÓT M, "AL Heé—i 5*1 —— 
for 4; M, Ex —— Lhar 1 E 1 (second interval) 
On, = =~ mi az, 1| 20 T) to (1 lt £M | are costo} — + 
d arc cos (M,— f£) + EEUU 
E 42—22 Mato QM 384-413) MA] 
(ale (105b) 
and tor 7 p- <t (third interval) 
0. -—5 (105c) 


Numerical results are presented in Part IV. 


INDICIAL LOADING FOR SINKING TRIANGULAR WING WITH SUBSONIC 
LEADING EDGES 


The six regions.—As in the study of supersonic-edged 
triangular wings, there are also in the case of triangular wings 
with subsonic leading edges various regions in which the ana- 
Figure 46 
shows the regions into which the subsonic-edged triangular 
wing can be most conveniently divided, where the trailing 
edges are again assumed to be supersonic and the solutions 
are cut off appropriately. Most of these regions have 


counterparts on the supersonic-edged wing shown in figure 37. 
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To begin with, region 6 lies within the spherical wave which The equations which relate the normal components to 
started’ at t=0 from the wing apex. Region 1 is within the | those in the free-stream direction have already been given 
cylindrical wave which was started at £—0 by one edge of the | in the section on region 4 of the supersonic-edged wing. _ 
wing, but outside the wave started by the other edge. Re- | Use of these relations leads to the following expression for — 
gionY£ is the area formed by the overlapping of thé two | loading in region 1 (in the coordinate system of fig. 47): ` 
cylindrical waves from the opposite edges, but outside the UT 





region influenced by the reflection of one of these waves on Ap\ Ba m Ms 41-4 m*rly|l—mz 

the opposite edge (secondary wave fronts shown in the P) =e mM tylim? “mMt—lylitme " 
figure). Region 5 is the area between regions 4 and 6 where TARNEN 

the flow is influenced by secondary (and higher order) ware are tan / pa tmu (106) 
reflections. Finally, regions 2 and 3 are similar to regions 2 Vivi-+m?+lyl—mz 


and 3 in the supersonic-edged case; region 2 being that unin- 
fluenced by. the starting phenomena and therefore having a 







loading already at its steady-state value, and region 3 being —7 \A erae M Dj l 

that which is unaffected by the disturbances emanating from dic NEU Do 0ynm(xaM, t) - DS 

the edges. — - y*-m(x +M,#) - ~~ : | 
i orooy extnh Mt : 


AN m*cinA 
.-- —8Starting spherical | fn fd -i | i 


L —— Leading edge now : Joan om yim 
——— Leading edge at t-a P os 
IIR ,.c-Secondary wave 
I] ue LAM reflections 


i .J7y7-mr- tym? r-y7mzr-t ylim? 





! \ ` x 
I \ FiGCRE 47.—Equations of Ines used in the analysis of the tMangular wing with 
k i subsonic leading edges. 
] a a ! M | 
Um PM Region 2: The loading on region 2, being the steady-state | 
^Slarting cylindrical wave loading on a triangular wing with subsonic edges, is well 


from right edge known. ‘The solution for region 2 of figure 46 is therefore | 


Fr;CRE £5,— The six regions used in the analysis of the triangular wing with subsonic given by (see, for convenience, reference 22) 


leading edges. -———— 
- . Ap eae 4a m* (x+ Mt) (107) 
Region 1: The solution for the load distribution in region 1 .N0Qo/: Eqm*(r4-M,— ^ 


Is the same as that for & two-dimensional wing starting | . | 
suddenly from rest and moving with a steady subsonic | where E is the complete elliptie integral of the second kind 
velocity V normal to the leading edge. A solution to the | with modulus J/1— 8?mt. EE 


latter problem for the initial part of the motion is presented Region 3: The loading in region 3 follows from reference 21 
as equation (A6(b)) m appendix A. In terms of the normal | and is 


components of velocity and distance, therefore, the loading (22 m. (109) - 
coefficient for the right-hand side of figure 46 can be written Qo/3 Ado 
immediately: 


Region 4: The loading in region 4 of figure 46 is calculated 


Ap, BU M, / i—zr, ae | by superposition, just as the solution for region 5 of the - 
C) eee ALN LFM, VM Mttr, 4 tan NV t—z wing with supersonic edges was obtained. The solution _ 
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in region 4 is tlie sum of the solutions for the right and left 
halves of region 1, minus the result for region 3. Thus 


mM, 


(22 Al ( tylim —y— me 
4 =, mM,tVi+m? v m Mist-4-y - mz 


"re pee MN ] dtum u 
Jim ) tare tan tvi--m VON mat 


r 


14 mi--y—mz 2 


Regions 5 and 6: In these regions the exact solution for the 
loading has not been determined. As was shown in the 
section on homogeneous boundary-value problems, such 
solutions would require the solution of a three-dimensional 
elliptic-type partial differential equation. <A later section 
will contain an approximate solution for these regions when 
the wing plan form is slender. 

Discussion.—AÀn isometric drawing of the load distribution, 
for the regions in which it is known, is shown in figure 48. 
Comparing the results for the loading on this wing to the one 
with supersonic edges (fig. 44), it is apparent that the prin- 
cipal difference in the two distributions is in the behavior 
around the leading edges; the loading being finite at the 
supersonic edge, whereas it becomes infinite at the subsonic 
edge. In view of the known steady-state results this differ- 
ence was to be expected. Elsewhere the loadings are quite 
similar. 


ERN MES Traces of regions shown in 
figure 46 
(Faramefers chosen so region 
5 is nonexistent) 





FIGURE 48,—Dfístribution of indielal loading on right panel of subsonic-edged tri- 
angular wing 


(25) m. Se (mM, „JEE 


The results presented in equations (106) through (109) will 
next be examined in a different light. Choose a given span- 


wise section on the wing and watch this section as time pro- 


gresses from (—0. This amounts to fixing the axis on the 
body and can be accomplished simply by using the quantity 
8 introduced in equation (93), 


s=m(z+ My). 


It is dar that 8 is the semispan of a given spanwise sections 
and that if equations (106) through (109) are written in term, 


of s, y, and i, for a fixed s they represent the variation of 


loading on & given section as time progresses. 
If the notation is further simplified by mtroducing the 
parameter 8, where 


1 
= - c. 7 (110 
equations (106) through (109) can be written in the following 


ad 


|. s—iyl ) 
are tan (8, —e4- yl (111) 
Ap 4ams 112 
Jo D 
(G2) =; A. | (113) 


(22 4 ar At = 38 (n MB. Ts a 


- = Lm 


m My, NDS —C-rFarc tan Vari 


(114) 


is ais =i) 


The load distribution across any section is given by equa- 
tions (111), (113), and (114) from the time ¢/8,=0 to (é/8,),, 
where the term (é/8,), is equal to 2s or s/m(44,+1)8,., which- 
ever is smaller. (At t/8,—2& the secondary waves shown in 
fig. 46 have just reached the spanwise section, and at (/8,— 
s/mB,(A4,+1) the spherical wave which started from the 
apex has just reached the spanwise section.) From (é/8,); 
to (£/B.)a=s/mB.(AL,—1), the loading has not been determined 
and from i/8,— (t/8))s to i= œ the loading is the steady-state 
value given by equation (112). Figure 49 shows this initial 
and final load variation plotted as a function of the parameter 
i/8,. At the beginning of the motion the loading is constant 
across the span, but this type of distribution is quickly 
modified arid the shape of the curve tends toward the steady- 
state loading given by equation (112) and shown in the figure 
as the distribution at é/8,=(t/B,)2. In fact, if the span is 
crossed first by the secondary waves rather than the spherical 
wave, when this span has traveled a distance such that 


-- M 
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Fite 49,—Vartation of indicist loading with 4/8, at a fixed spanwise section. 


i/8.=2s, the expression for the loading given by equation 
(114) becomes 

(= 
. Fo 


_ 16m. se 


n » {11 
er (115) 


— 
which differs from the value given by equation (112) only by 
a constant of proportionality. Both before and after the 
time £/8,—2a the shape of the loading curve varies from ia 
simple type represented by equation (115), but the trend is 
established. 

The average chord loading factor P$, introduced by equa- 
tion (92a), can now be determined. for certain regions. 
Hence, if the notation 


—(tl8 
-x Sp dy 
P= Sza —s Qo 


is adopted, there results for the early part of the motion, 
that is, for O< r/B. (7/84): 


(116) 


p,-2 (2 -Z)+ 4mB,Mo(r/B,) (117) 


Equation (117) was derived by integrating equations (111), 
(113), and (114). For values of 7/B,2— I1/B,m(M ,— 1)=( 5) 
equation (112) is valid. Hence for — aoe 8 ) 


[d 


P 2r m Mo 


CEP (118) 
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Figure 50 indicates the magnitude of this average load for 
both large and small values of 7/B,. Notice that for small 
values of r/B, it is sufficient for the establishment of the 


curve to specify the parameter mM,8,, but for large values - 


an additional parameter must be given (such as Afa in the 
figure). Notice, further, that in spite of the large variation 
in the distribution of the loading, as shown in the previous 
sketch, the average value P, varies linearly throughout the 
intervals considered. This result is similar to the one 
obtained for triangular wings with supersonic edges. 


4 
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FIGURE 50,—Variation of P, with (Be at a fixed spanwise section. 


INDICIAL LOADING ON YERY SLENDER TRIANGULAR WINGS 


In Part II of this report it was pointed out that if the wing 


is slender (i. e., has a small ratio of span to chordwise length) 
the governing partial differential equation simplified to the 
form 





Pu— Py — P= 0 (119) 
where the independent variables refer to a coordinate system 
that is either fixed on the wing or is fixed with respect to the 
still air at infinity. The boundary conditions that apply 


when the axes are fixed on the wing will now be considered _ 


in some detail. (In Part II the variables z;, yi, z; and i 


were used-to denote this coordinate system. In order to 


avoid a cumbersome notation, however, the subscript 1 will 
be deleted in the following. ) 

Just as in the previous sections of this report, consider & 
triangular wing which is at rest for <0, starts suddenly to 
move at a forward velocity equal to V, at (—0, and continues 
at this same velocity for £750. In this case, the triangular 
wing can be considered as having & finite chord, since span- 


wise sections act independently in slender wing theory, and 


the final integration in the x direction can be stopped at any 
desired chord length. It should be emphasized that in this 
case, To may be either subsonic or supersonic. A section in 
the spanwise direction, as for instance section AA in figure 
51 has & trace in the yé plane which is a narrow rectangular 
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FIGURE 61.-—Blender wing In ry plane and trace of section In pt plane. 


strip along the taxis. Since equation (119) has been derived 
on the assumption that the velocity gradients in the y, z, 
and £ directions are independent of the gradient in the z 
direction, the boundary conditions along the strip shown in 
the figure are independent of those on other strips corre- 
sponding to traces from spanwise sections along the wing. 
Hence, the problem is to find a solution to equation (119) 
which will make e, constant over the strip (since a given 
spanwise section experiences constant downwash whether the 
slender triangular wing is undergoing sinking or pitching 
motion) and at the same time will satisfy the other condi- 
tions listed under equation (2). In the lifting-surface 
analog this corresponds to the problem of finding the velocity 
potential over a flat rectangular wing of low aspect ratio 
situated in a free stream moving at à Mach number equal to 
42. Solutions to the latter problem can be obtained by 
various techniques, and so the procedure will be first, to find 
the potential for the steady-state, flat, rectangular wing, arid 
then, bv analogy, to convert this to the solution for either the 
sinking or pitching slender triangular wing. 

The steady-state, lifting-surface problem.—Lift.ng-surface 
solutions for the loading on & rectangular wing traveling at 
supersonic speeds have been developed for regions 1, 2, and 
3 of figure 52 (by Busemann and others), and by means of 
these solutions the load distribution on a spanwise section of 
the triangular wing can be determined to a time necessary 
for sound to travel that span length. For £752s, however, 
the solution becomes considerably complicated by the in- 
creasing number of reflections from the edges. Reference 24 
gives solutions for the loading on & rectangular wing in 
region 4 and indicates methods for extending the solution to 


REPORT 1077—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


2s 





FIGURE 52.—Reglons used in the discussion of the low-aspect-ratio rectangular wing. 


regions farther along the wing. Already in region 4, how- 


ever, the expression is cumbersome and in higher-numbered 


regions the expressions become difficult to manipulate. 
These methods, therefore, will be discarded in favor of a 


more approximate but simpler analysis. 


If z is the distance along the chord, y the distance along 
the span, and ¢ the semispan, then the solution for regions. 
1, 2, and 3 of figure 52 ean be written (for convenience, see 
rarene 22): 


Region 1 
Ap w 
d 4 Y (12Ua) 
Region 2 
| Ap  8w SS | 2 
| o V arc tan ior (2 Ob) 
Region 3 
Ap 8w Ty —y T7 
dr arc tanq/ — — tare tan "gd rae > 


(120c) 


n o= Ae Loi. + ABuaL- 
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As x increases (1. e., for higher-numbered regions in fig. 52) 
it Is reasonable to assume that the spanwise variation of 
loading is virtually invariant with xr, that it is "smooth," 
and that it falls to zero at the side edges. 
fore, that the loading is given by the relation 


wy) aan 


Spanwise, this has the variation shown in figure 53; chord- - 


wise it is as yet arbitrary. To fix the chordwise distribution 
the value of f(r/s) will be determined so that the vertical 
induced relocity along the center line is constant and 
equal to w. | 

The solution to this somewhat artificial problem approaches 
the exact solution to the steady-state lifting-surface problem 
for a flat rectangular wing along sections far behind the lead- 
ing edge; closer to the leading edge it only approximates the 
exact solution; and, of course, in the vicinity of the leading 
edge it will be least representative. But, on the other hand, 
the exact solution is known in the vicinity of the leading 
edge and it turns out that the solution of the problem posed 
above forms a reasonable continuation over the remainder 
of the wing. 

The velocity potential for the problem which has been set 
ean be readily expressed in terms of an integration of ele- 
mentary horseshoe cue over the plan form. Since the 
Mach number equals 42, then according to reference 19, 


eS oe — 21) (áp/qgdzid ys 
(qt z2)4 (c — n — yi — 2? 





where A is the area on the wing within the forecone from the 
point P(z,y,z), at which ¢ is to be determined (the shaded 
area in fig. 53). 

The simplification of the last expression is given in refer- 
ence 19. The result is the integral equation 


L=fn)t= [Fada ndo. (122) 


where ?—z/e and G is given bv 


Ei L<a-m 
E4—(1— kK: 


L2 


G(»—7n)— 
7—7m $1 


l 
kı =——; ka= — 71 


7—"71 

The modulus of £i is $, and the modulus of K; and F; is ke. 

The solution of equation (122) for f(y) is not difficult 
when numerical methods are used. For intervals of s; equal 
to 0.2, the result is given in tabular form.in table IIT, and 
also in figure 54. As mentioned previously, the function 
f(g) determined by this approximate method will be least 
representative of the exact solution in the region near the 


Assume, there- ' 


Ap 








P(x,y,z) 





FIGURE 53.—Span wise loading and region of Integration used in the analysis of the low-aspect- 
ratis rectangular wing. 


leading edge. However, since the exact solution is known 
for x/s—=n<2, a compromise can be effected. If the span- 
wise average of the loading is calculated from equations 
(120), it is' found to be & linear function of «, starting at 
4/x for n=0 and falling to zero at ņ=2. On the other hand, 
the spanwise average of loading given by equation (121) is 


— Y $n). Therefore, an improved solution for f(y) consists 
0 


in taking /(q) as the value given by 


2 

f= 2-2) 

in the region O<y<2, and then fairing this curve into 
that given by the solution of equation (122) for y>2. The 
two curves are shown in figure 54. 

By using the results listed in table III, the loading over a 
low-aspect-ratio rectangular wing flying at a Mach number 
equal to 4/2 can be estimated. Of particular interest is the 
dam ped oscillatory nature pf the load, as shown in figure 54, 
falling to zero &t one span length behind the leading edge 
and taking alternately negative and positive values beyond . 
this point. (See reference 25.) A somewhat different 
approach to this problem (reference 26) has led to a solution 
very [ike the one given here. 
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TABLE HIL—SOLUTION OF EQUATION (122) necessary, therefore, to operate further on the solution given 
p | “|..| for the loading in the steady-state problem to obtain the 
= | solution for the loading in the unsteady prout But 





`- "t" mA 











p T [ ate] = xl ici 

4 44 

6 46 

.8 48 "t NOE 
ri 52 so that if the notation 

4 bc M Li 
(OH T £m ML 
1 20 6.0 p - 
| 2.2 8.2 e 

24. 8.4 E f 
i x zi is adopted (where a represents the loading in the anal- 
| i T$ ogous steady-state problem), then the expression for the 
Lus ced A 

unsteady loading can be given in terms of (22) by the equa- 
$ 07s 

fom tion 


rw nit P) tmo: | (= ) | (124) 


„From equation (120) 








-From equation (122) By the ney equation (124) to equations (120), 

d the loading for the various regions in the yt plane of the un- 
steady wing can be found. For region 2 in figure 55 there 
results 

Em °( p b De 
age A arc tan AL MAN Hagan 
r 8a, are tan EUN dt jl 
FrongE 54.— Variation of fín) with s. - t- |x| T É —8 8-- iyi 


E ; : which becomes 
The unsteady-analog, sinking wing.— The first step in Nm 


deriving the unsteady-flow results for the sinking wing from $a i—e+ly u a—|y 

the steady solution is to replace x with £. In equation (122) ur M, |m M of z]y] Ls are tan y E. | 

this corresponds to replacing 7 with + where r is equal to - 

t/s (equations (116)). The second step is to replace w/V, The lodding coefficient can be similarly derived in the other 


with —« and to rederive the expression for loading coeffi- | regions so that finally, for the regions shown in figure 55, 
cient since In the time-varying problem it is expressed in a 

somewhat different manner than in the steady-state analog. p 

In the unsteady case, as the triangular wing moves through Section AA 


& fixed reference plane the local span intersecting this plane 
grows as a function of time and. equation (12), which repre- 
sents the partial derivative with respect to time with x fixed, 
must be expanded to the form 


wovon otk van at LL os dac) 
Go n Voll, ot VAL, ot 


where | Se and [55 indicate -derivatives taken at 
1 t 








x 
constant s and É, respectively. | Since 8 is n to m(x-- Aa t}, Figure 55.—Reglons used to express loading on section of unsteady, subsonic-odged, tri 
Os/Of equals mAf,, and there Restle | EN , angular wing. 

<P. (861 mM Ie] ) Region 1 
123 
A, eae um "D ^p 4^ . (95) 


In the steady-state problems an analog to the term | Region 2 


involving El is missing, and the loading coefficient is Ap Aen M [ize He -Farc an IL) mE 
Go do M, : 7 i—s+iy| 


given entirely by an operation sautvdsnt to [Ax]. It is (125b) 
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Region 3 
Pa Sa Eo i—3—y 
Se (m mM s ema i sy 5 
are tan [Ever tan 1 [etx 0) 


(125c) 


For the interval £22s equation (121) must be considered. 
By means of equation (124), the expression for the loading 
coefficient can be written 


Pay. M, E Cw A ) ena f. ta t oN (Hy at | 


whieh becomes 





Ap m M, 


Pur (1—mMor)4/ 1— (¥Y fo + 75 
V, il = Ju) 





1i ye 


(126) 


where f(r) is the solution to equation (122). Notice that for 
large r (when the loading has reached its steady state), f(r) 


is zero and [ j f(rüdriisunity. (See table III.) Hence, the 
loading is given by the equation 


Ap gams 
do 482—347 





which is the steady-state value for a slender triangular wing 
(equation (112) when £-—1). 

It is now possible to derive the chord loading factor P, 
as defined bv equation (92a), thus 


P _ Mo L AP dy 


T 
2sa —t Qo 


Placing equations (125) and (126) in this expression, it is 
found that for 0<7r<2 


Eu 2(2— r)d- £m Uar (1218) 


and for 722 
P,—s1—mMsr)fir)-2mM [' fe)dr (127b) 
af) 


Since the values of P, given by equations (127a) and (127b) 
were derived using different methods, their magnitudes at 
r—2 are not equal. The final curve for P, must be con- 
structed by fairing the solution for 7 €2 into that for 7772. 
Figure 56 shows these results together with the final curve 
chosen (solid line). 

The unsteady-analog, pitching wing.— When the wing is 
pitching at a steady rate about its apex, the equation for 
the vertical induced velocity on the plan form is 


p,y— —(z-.- Mg 
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ov 
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c —£fquation (I27& 
» jd /e7 5 





G 2 4 6 8 j m. 
r Es 
FiGCRE 56.— Varlation of Pg with r. 


so that —w/V, in the steady-state equations (120) and (121) 


becomes ĝe/mV. Since the loading coefficient is still given 


by equation (124), there results for the conversion of equa- 
tion (121) the expression | 


oo re | ey 
Qo ~ MAL ao M 


and this can be reduced to the form 


yf) mM ` Tna N "Kadr | 


Ap__ 46 — 
a RT A EE 


LAE ta i) L Weide 


As in the discussion of equations (127), it can be seen that 
equation (128) becomes for the steady state (7 large) 


mr) s? ey fi) 





(128) 


Ap 40 (257 y 

do Voi x's? y? 
and this can be shown to agree with the steady-state slender- . 
wing results given in reference 27. 


It is now possible to derive the chord loading factor P, 
as defined by. equation (92b) 


| m.MaVo r ÂP dy 
2g? -t mt 








Using equation (128), one finds for 7752 


P,—a«(1— mM) f(r) +3m Mor i f(r)dr — (1298) - 


and & similar analysis based on equation (120) vields for 
0€7€2 u 


P zt (2 — T +4mAfor -3 mMor*) (129b) 
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As in the case for P4, the two equations for P, do not join at 
r—2 and the final curve must be constructed by fairing the 
soluuon for 7 €2 into that for 722 (see fig. 57). 






l r--Equation (129b) 
P, 2 2^ , Equation (/298) 


f f 
fyé 


FIGURE 5:.— Variation of P4 with r. 


Discussion.-—It is now possible to assess the accuracy of 
the solution for verv slender wings in the interval 0€ 7 €2 


by comparing the values of loading given by slender wing - 


theory (equations (125)) with the exact loadings given by 
equations (111), (113), and (114), and also by comparing 
the slender wing value of the chord loading factor Pj (equa- 
tion (127 a)) with the exact value given by equation (117). 
It is apparent that the approximate solution differs from the 
exact only by a stretching factor in the ¢ direction. Hence, 
if r is replaced by 7/8, and m (note m is proportional to 
y/t) by MBa where 8, is given by equation (110), then equa- 
tions (1258), (b), and (c) are identical with equations (113), 
(111), and (114), respectively, and, of course, equation 
(127a) corresponds to equation (117). 

'This rather remarkable result can be enlarged upon from 
another viewpoint. Suppose that in the steady-state analog 
problem the wing had been flying. at some Mach number 
other than 4/2, say Ae The solution to such a new problem 
could be obtained from the old one merely by applying the 
Prandtl-Glauert correction, that is, by stretching all dis- 
tances in the x direction by the factor 1/8, where 8?=|1— 
Mèl. Such a procedure would convert, for example, equa- 
tion (1278) to the form 


Po8,—2(2— rfp) +4m B.M, 7/Be 
Finally, if P, is adjusted so that P;—4 at 7—0, there results 
P,—2(2— r[B,) 3- 4$ Msc 


which is exactly the answer given by equation (117). It is 
possible to simplify the statement of the procedure by simply 
remarking: The exact results for Ap/qo or P, in the interval 
0Xr7/8,S 2 can be obtained from the approximate results 
for a very slender wing by making an effective Mach number 
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correction to the right-hand side of equation (125) or (1272), 
respectively. 

It is interesting to pursue this concept even further. 
Consider a spanwise section of a triangular wing as time 
increases from the starting impulse. The primary wave 
fronts emanating from either side pass across the section, 
forming the Mach lines in the steady-state rectangular-wing 
analogy. For very slender wings these lines make a 45° .. 
angle with the trace of the side edge and are used to divide 
the plan form into regions as in figure 58.. Now find the. 
actual position of these primary wave fronts as they form a 
trace on the section in the yr plane. A straightforward 
calculation shows that these lines actually make an angle 
equal to arc tan 1/8, with the trace of the side edges. Hence 
the effective Mach number which is used to correct the 
slender-wing results in the interval 0& 728, is that which 


" makes the Mach lines of the steady-state analogy coincide 


with the actual trace of the primary wave fronts. 


r--Acfual trace of 
primary wave front 


-m ee oe oe 






Trace of primory wave 
front obtained from . 
slender-wing theory -- 


FiaUuRE 58,— Wave traces on spanwise section of slender wing. 


` INDICIAL LIFT AND PITCHING MOMENT ON VERY SLENDER TRIANGULAR 


The lift coefficient for the sinking wing is given in the 
notation introduced in equation (92a) by the equation 


co 2 
Cas f 8 Pdr 


mejo Mo 


where P, has been determined in the last section as a function 
of r=t/s and ¢ is the root chord. Consider the situation 


at a certain fixed time and let the x coordinate in the above _. 


formula be fixedin the wing. Then set 


=- c (180) 
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and as before 


_M 


ro 5m Matr (131) 


where r, is the number of wing-chord lengths traveled. In 
this way it is possible to obtain the expression 


1 
AVAGA =2 Í tP, =H) dt + (182) 
and similarly 
M,C =-2 f 2P TAL (133) 
irt US ma T ix 0 m Mat ie 


The equation for the lift and pitching-moment responses 
(where again the pitching moment is taken about the apex) 
on & pitching wing are s 


MC: 


t int : To ' 
M, "(s =2 | eP (zg) 039 
T, 
and 
1 f I 
MC, (= CR 2 DP, zu) dt (135) 


Eo 


0 


The values of Pa and P, were taken from curves similar to 
figures 56 and 57 in the last section (using the faired curves 
in the vicinity of 7—2) and the results for the indicial lift 
and pitehing moment in terms of fẹ the number of chord 
lengths traveled, are shown in Part IV. f 


PART IY—RESULTS AND DISCUSSION 


TWO-DIMENSIONAL RESULTS 


The methods presented in the previous sections have been 


used to calculate the indicial lift and moment curves for two- 
dimensional wings flying at Mach numbers equal to 0, 
0.5, 0.8, 1.0, 1.2, and 2.0. These results are presented in 
figure 59. | ] 
Figures 59 (a) and (b) show the variation of the indicial 
lift and moment on a sinking wing. "The initial value of the 
lift is given bv the expression 4£o/M,. When the free- 
stream Mach number goes to zero this expression still is 
valid, the initial value being a pulse of force that occurs at 
£—0. . The final values of the lift are simply the two-dimen- 
sional, steady-state results given by the Prandtl-Glauert 
rule. Figures 59 (c) and (d) show the lift and moment 
yariation on @ wing pitching about its leading edge. These 
functions were not computed for 14,—0.5. The results are, 
of eourse, subject to the restrictions of linearized, com- 
pressible-low theory and, for example, the calculated 
responses given in figure 59 for sonic speeds must be con- 
sidered as being outside the realm of validity within a few 
chord lengths of travel. In application to high-frequency 
oscillations, however, the initial portions of the indicial 
eurves dominate the response characteristics of the airfoil 
and calculations near M, equal to one need not be invalid. 


More accurate values of the indicial curves are given in 
tables I and 1I for the subsonic Mach numbers 0.5 and 0.8. 
These tables can be extended to larger values of rọ by means 
of equations (71) and (72). . 





e 4 6 
Chord lengths traveled, V,£'/ 


(a) Lift on a sinking wing. u 
(b) Pitching moment (about leading edge) on a sinking wing. ` 
FiGCRE 59.— Variation of two-dimensional Indfcial lift and pitehing-moment coefficients with 
chord lengths traveled forseverai Mach numbers. 
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dimensional wing only at the beginning of the motion and 
again when the steady state has been attained. 

‘Second, notice that, since all the characteristics. for the 
triangular wings are independent of the angle of sweep, they 
are vàlid for any unyawed triangular wing as long as the - 
edges are supersonic. 

Third, it is apparent that the transition of the total indicial 
lift from. its initial to its final value is less abrupt for the 
triangular than for the two-dimensional wing. The initial 
.&nd final values of C; depend on 1/Af, and 1/8, respectively, 
so that as the Mach "number i is increased the variation dies 
out altogether since 8 and Ma approach one another. The 
same remark applies to all the other coefficients. 


m 


MARO S D IER TRIANGULAR WINGS EM ES 

The results for the lift and moment on a slender triangular 
wing that is sinking or pitching about its apex are shown in 
figure 62, The analysis by which the results were obtained 
is valid when both m, the tangent of the semiapex angle, and 
mM, are small; the results are given for mAM,—1/8. The 
_ curves are all qualitativ ely alike, in each case the response 

falls from its high initial value to a minimum at about rg— 1/8 | 
and then recov ers and practically attains its asymptotic 


value at =l. 
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Chord lengths traveled, V t'fe 


(c) Lift on wing pitching about its leading edge. 
(d) Pitching moment (about leading edge) on a wing pitching about fta leading edge. 


FiauRE 59.—Coneluded, 


TRIANGULAR WINGS WITH SUPERSONIC EDGES 


Curves of the indicial lift and moment on sinking and 
pitching triangular wings with all supersonic edges are shown 
in figures 60 and 61 for Mach numbers 1.2 and 2.0, respec- 
tively. For the purpose of comparison, the curves for a 
two-dimensional wing flying at the same Mach number 


have been included in each figure, as well as the curves for | ELT. i | 

the same triangular wing in reversed flow. C 
Several conclusions can be drawn from these results. ‘ol io] |. 2 ae: | 
First, notice that the total. indicial lift on the triangular ; NEN" 
e 2 6 8 


sinking wing is the same at every instant as that on the same 4 
od t Chord lengths traveled, Vat/e 
wing in reversed flow (both wings, of course, having started. 


with the same velocity at the same time), and that the value (a) Lift on sinking wings. 


a a a . ‘oe Li F G .— Indi ` s h 
of this lift is the same-as the total indicial lift on the two- Bc acca ka aa a ae 
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«£f E 
Chord lengths traveled, V,t'/c 


ib) Pitching moment (about leading edge or aper) ou sinking wins. . 
Figure (0.— Continued, 
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Chard lengths traveled, Vo t'/e 


(c) Lift on wings pitching about their leading edge or apex. (d) Pitching moment (about leading edge or apex) on wings pitching 
about their leading edge or apex. 
Figure 60,—Concluded. 
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“O 4 .& Leg L6 Z0 >- 24 Q 4 8a 3 Ó3ómaAe L8 2.4 
, Chord lengths traveled, V, tjc Chord lengths traveled V, tfc 
ia) Lift on sinking wings. (c) Lift on wings pitching about their leading edge or apex. 
(tb) Pitching moment (about leading edge or apex) on sinking wings. (d) Pitching moment (about leading edge or apex} on wings pitching about their leading 
FIGURE 6t.—Jifdiclal aerodynamic characteristics of triangular wings with supersonic edges; edge or aper. 


Mis 2.0, FIGERE 01.—Concluded. 
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FIGURE 62,—Indiclal aerodynamic characteristics of smking and pitching triangular wings with slender plan forms. 
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2 4 o) 
Chord lengths traveled, V tc 


ta) Lift, Mem 1.2. 
(hb) Pitching moment (about leading edge), Mim 1.2. 
FIGURE 63.—Indiela] lift and moment on sinking rectangular wings fiying st supersonic 
epeeds. 


RECTANGULAR WINGS 


By means of the analysis and curves presented by J. W. 
Miles in reference 2, the indicial results for a two-dimensional, 
sinking wing can be compared with those for a rectangular, 
sinking wing flying at supersonic speeds. This comparison 
is given in figure 63 for Mach numbers equal to 1.2 and 2.0 
and aspect ratios equal to 2, 4, 6, and infinity. It can be 
seen that the modification in the shape of the indicial curves 


212183—54———29 
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& 
Chard lengths traveled, V, tfe 





-g a g 
Chord lengths traveled, V; tc 
(c) Lift, X422.0. 
(d) Pitching moment (about lesding edge), Xe. 0. 
FiovRE 63.—Concluded. 


Le L6 


brought about by the decrease in aspect ratio is similar to 
that caused by a decrease in Mach number from supersonic 
to subsonic magnitudes. 


AMES AERONAUTICAL LABORATORY 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
NlorreT? FELD, Cauir., October 12, 1950. 
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APPENDIX A 


DETERMINATION OF SUBSONIC, INDICIAL, 
SECTION LIFT AND PITCHING-MOMENT CURVES 


THE LOAD DISTRIBUTION 


The following results for the indicial load distribution on 
sinking or pitching wings can be obtained in two ways. 
One of these methods will be outlined in the subsequent 
paragraphs. The other is outlined in references 28 and 29 
and is referred to as the lift-cancellation technique. The 
latter method has been used to check the load distributions 
originally obtained by the former so that an independent 
check of these results has been carried out. 

It was shown in Part II, Two-Dimensional Boundary- 
Value Problems, that the lifting-surface analog to the solu- 
tion for load distribution over an unsteady, two-dimensional 
wing traveling at a constant subsonic speed involved the 
calculation of load distribution over a swept-forward wing 
tip with subsonic edges. Figure Al indicates the geometry 
associated with the boundary conditions. Solutions are 
given only for the loading in the five regions shown. As has 
been discussed in the text, for Ma=0.8 this was considered 
adequate to define the behavior of the indicial responses in 
the early stage of motion. For M,—0.5 the loading in regions 
farther down the wing had to be calculated; and in all, for 
this Mach number, the loading was analyzed in the 11 
regions shown in figure 25. However, the latter analysis 
was carried out by a numerical application of the lift- 
cancellation technique and none of the details will be pre- 
sented here. 

In the notation of the unsteady problem the expression 
for the velocity potential can be written 


-—Í hi K m Ch a zo 


where r is the area on the wing plan form included in the 
Mach forecone from the point (¢,7). Equation (Al) is 
applicable only for cases in which w, is known at all points 
within the forecone, as is the case when the edge of the wing 
within the forecone is everywhere supersonic (1. e., region 1 
in fig. Al). However, Evvard (reference 30) has. extended 


-(A1) 


the solution provided by equation (A1) to include cases such 


— 





Figure À1.— Regions used in analysis of unsteady, two-dimensional wing flying at subsonic 
speed. 


as shown in figures A2 and A3 in which the forecone inter- 
sects a subsonic edge and includes a region of unknown 
upwash. As was pointed out in reference 30, equation (À1) 
applies in these instances if the area of integration r is limited 
to the shaded regions shown in the figures. It is apparent, 
therefore, that the potential (and thus the loading) over a 
sinking or pitching wing can readily be determined for regions 
1, 2, and 4, in figure Al. 

Points in regions 3 and.5 have forecones which intersect 
two subsonic edges, and the method just discussed can no 
longer be directly applied. In reference 30, however, a 
method was given of evaluating the upwash in the region 
between the Mach cone from the apex and the leading edge 
(region 6 in fig. A4). Thus, the plan form has become, 





FIGURE À2.—AÀres of integration when Mach forecone intersects subsonic trailing edge, 





t 
FiouRE A3,—Area of integration when Mach forecone intersects subsonic leading edge. 
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FiaURE 44,—Region 6, area between wing leading edge and trace of foremost Mach oone, 
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effectively, one such as shown in figure A5 in which only one 
edge is subsonic. This reduces the problem of finding the 
potential in these regions to the same problem as was involved 
in region 4. The analysis used in finding the loading over the 
various regions will now be considered. 
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FiGrR£ ÀA6.—Reglons of figure A1 in the rs plane. 
First, introduce a new coordinate system in which the lines 
z- —i and r=i are taken as the r and s axes, respectively. 
(See fig. A6). This amounts to s rotation of the original 


system of axes through an angle of 45°. The transfor- 
mations relating the r,s to the x,y system are 


r= —r) pnm 4 (r48) 
l 
=; (t+ x) = (s—r) 
In the new coordinate system equation (A1) is written 
Wari, 8)dridai &,)dr.ds, 
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The vertical induced velocity w, over the wing plan form 
is given in equations (2) and (3) of the text for the sinking 
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wing and pitching wing, respectively. The method developed 
in reference 30 was used to obtain the value of w, over the 
area between the lines s=0 and r—s(1-4-MQ/(1 —Mq) (region 


6 in fig. A4). The results for the sinking and pitching wing 
are respectively, 
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In the r,s coordinate system, if the subsonic trailing edge is 
not included in the forecone from the point at which the loading 
is being determined, the expression for the loading coefficient 
can be written 


io cT. (orae) | ka 


However, if the trailing edge is included in the forecone, 
and if the Kutta condition is to apply along such an edge, it 
can be shown that the equation for the loading coefficient 
assumes the form 
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The preceding method can be applied to the sinking wing 
to obtain the following integral relationship for the loading 


=- - * A * >. 
over the various regions. The subscripts on (25) indicate 
0 
the region for which the particular equation applies. 
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TWO- AND THREE-DIMENSIONAL UNSTEADY LIFT PROBLEMS IN HIGH-SPEED FLIGHT 


LIFT AND PITCHING-MOMENT COEFFICIENTS 
The lift and pitching-moment coefficients may be obtained 
by suitable integrations of equations (A6) and (A7) and 
are given in the time intervals indicated in figure 34 by the 
following expressions: 
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Pitching wing 
SS TM, LE 


etr p e 5-3 My F&041)41 | (A108) 
f. Mg 


r Saf lal i 3,201 arya 
ges AL Ee MJ 5 My Hita My 


pal) +1} (A118) 


1 1 
1 TAS SET — M, 
3) 4 


a fo(—5 Me? +24f,— 
ere T9] srra Md 
 O304My —\ == Tagg) + 

CAN SR i| 


(vd 
VETSMTA EC +5) 
to + I — tto I Me t h^! 
[es en Viri I 24-MPC | 
M 
(x a hf m wet: dz, (Al0b) 


El 4M, i+ 
where G is defined under equation (A7e). 
E. 8 mE e 


T = aM 341—Afà — Mj 
1—45 ys «1 HAL} — il- 
VANILA LAGO -MF 


2 
| "a; Oe DC Méca Mgr 


arc IER 








; anh" 








fos —ary—2 t MEN Mu 
git — M9 1 M SN &(1J-Mg—1 
Mg 1 t t(l — My +1 
s 63) arc an AETI SU MO I 


ù 





{t MYB MA | 609349—634,—11) | 
[240 FI) T20 FM 


3941+ 541f,+19 
SMS J 9-0 Ma 


[ 464345 tees VoM tau gn " fc A [E c 


(A11b) 





446 


REPORT 1077—-NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX B 


DETERMINATION OF SUPERSONIC, INDICIAL, 
SECTION LIFT AND PITCHING-MOMENT CURVES 


THE LOAD DISTRIBUTION 


In the case of the unsteady supersonic wing the expression 
for the velocity potential may be readily obtained by placing 
the values of w, given by equations (2) and (3) of the text 


in the equation 
=-if | aj Wull £ 
V~ — (z — zi Tı) 


where 7 is the area on the plan form included in the Mach 
forecone. The loading may then be calculated from the 
relationship given in equation (12). 
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SINKING WING 


The load distribution over the regions A, B, and C shown 
in figure 34 are given by the following expressions: 
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PITCHING WING 
For the case of the pitching wing the values of at 
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in regions A, B, and C are, respectively, 
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LIFT AND PITCHING-MOMENT COEFFICIENTS 


The lift and pitching-moment coefficients may be obtained 
by suitable integrations of equations (B2) and (B3) and 
are given in the time intervals noes: in figure 34 by the 
following expressions: 
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